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PART I. THE MOTION OF THE POLE 
NUTATION AS DERIVED FROM LATITUDE OBSERVATIONS 


By E. P. FEDOROV 
Gravimetrical Observatory of the Academy of Sciences of the Ukrainian S. S. R., Poltava, U. S. S. R. 


Abstract. The results of several long series of latitude observations have been used for a separate determination of the 
coefficients of some nutational terms in obliquity and longitude. 


The derived value of the constant of nutation is essentially smaller—and those of the coefficients of the semiannual and 
semimonthly terms larger—than the respective theoretical values, determined on the supposition that the earth as a whole 


is an elastic body. The theoretical value of the ratio of the axes of the nutational ellipse evidently does not need any 
correction 


On the basis of these data some conclusions of a qualitative character are made on the interaction between the core and 
the shell of the earth. 


The object of the present paper is to show that investigation of nutation. For this purpose we 
some conclusions about the interaction between must derive directly from observations more de- 
the earth’s core and shell can be drawn from an __ tailed information concerning nutation than was 
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hitherto available. Latitude observations over 
long intervals are likely to give the most favor- 
able material. The main, semiannual and fort- 
nightly terms of nutation deserve special at- 
tention. 

The effect of the main term on declinations 
may be expressed by the formula: 


Aod =— No (mo cosasin® — sinacosQ), (1) 
where 


No = the adopted value of the constant of 
nutation, 
ny = that of the ratio of the axes of the nuta- 
tional ellipse, 
a = the right ascension of a star, 
Q = the longitude of the ascending node of 
the moon’s orbit. 


Let us suppose that both No and 7» need cor- 
rections, and besides that there exists a lag of 
phase different for nutation in longitude and 
obliquity. Then equation (1) should be replaced 
by the following: 


Aé =— (Ny, + AN) 
X [ (mo + An) cos a sin (Q — B) 
— sina cos (Q — Be) ]. (2) 


If the actual effect of nutation is expressed by 
formula (2) the analysis of latitude observations 
should reveal the difference 


Ag = Aod — Ad (3) 


since in the reduction of observations formula (1) 
was employed. This difference may be repre- 
sented in the following form: 


Ag = Aicosacos2 + By, sin a cos 2 
+ A.cosasin2 + Besinasin®, (4) 


where 


Ay a Nonobx A» = NoAn — ANnpo ¥ 
By, SS AN Bs ae NB2 j (5) 


Thus, our problem is reduced to the determi- 
nation of the coefficients A1, A2, Bi, Bs. I decided 
to use for this purpose observations of the inter- 
national latitude stations. However, the results 
of these observations, as taken directly from 
publications of the Central Bureau, are unsuit- 
able for an analysis for deriving the periodical 
variation of the form (4). It is necessary first to 
apply certain corrections. This was done by T. 
Hattori when he used the latitude values derived 
at the international stations for the determina- 
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tion of the constant of nutation (Hattori 1947 
1951). My preliminary calculations in some re: 
spects resemble those-of Hattori but they differ 
substantially in some points. 

Earlier I noticed that sometimes the scak 
values adopted by the Central Bureau of the I.L.S 
had been subject to considerable errors. Con 
trary to Hattori, | made an attempt to free the 
observed latitudes from the effect of these errors 
as well as that of a variation of the mean lati 
tudes. Thanks to Uemae’s work (1953), it be 
came possible to exclude the errors made by the 
Central Bureau in applying corrections for th 
Ross terms of nutation. 

I do not dwelfhere on the description of eam 
phase of my calculation as all the necessary de 
tails are given in another paper (Fedorov 1958) 
I confine myself to giving the final result of thi 
calculation. From the analysis of about 135,00 
observations at Carloforte, Mizusawa and Ukial 
the following expression has been obtained: 


Ag =— 070081 cos a cos 2 | 
=e25 | 

— 0”0128 sinacos2 + o*0120cosasinQ | 
£25 +19 | 

— o%0004 sinasinQ. (6 

+19 | 


The same material has been used for derivin, 
the fortnightly term in latitude variation, but i1 
addition I have availed myself of the result ob 
tained by H. R. Morgan (1952) from observa 
tions with the Washington P.Z.T. from 1931 t 
1951 and that obtained by A. J. Orlov (1952 
from observations with the Pulkovo zenith 
telescope from 1915 to 1928. Thus, the tota 
number of observations used for deriving thi 
fortnightly term is 230,000. Some details of thi 
calculation are given in two other papers (Fee 
rov 1955, 1958). The final result is: 


Ag =+ 070086 sin (2 © — a) |) 
+14 a 
— o%0019 cos (2 € — a) 
+6 
+ o%0021 sin (2 € + = 
=) 


+ o”oo01 cos (2 € ae (7 
+10 


It will be of much interest to have also ai 
expression for the semiannual term. However, it 


~ this case the analysis of routine latitude observa 


tion meets special difficulties. The only observa 


ie 
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tions which are likely to provide favorable data 
for such an analysis are those of two bright zenith 
stars at Poltava. An attempt to derive the semi- 
annual nutational term from these observations 
was made by N. A. Popov. He obtained the fol- 
lowing result: 


Ag = 07027 sin (20 — a), (8) 
+4 


© being the mean longitude of the sun. 

Using the results of (6), (7) and (8), I have 
derived the expressions of the three above-men- 
tioned nutational terms both in longitude (2) 
and obliquity (€s — «o). They are given in the 
Table together with the theoretical expressions 
of the same terms. The latter have been obtained 
assuming for the constant of nutation the value 
9"220 which was found from the following rela- 
tion between H, the mechanical ellipticity of the 
earth, uw the ratio of the moon’s mass to that of 
the earth, and the constant of nutation: 

N = 231981"8H ees he (9) 

I+p 

Both this relation and the theoretical expressions 
for the nutational terms given in the Table were 
first deduced on the assumption that the earth is 
rigid, but they would be practically unaffected if 
allowance were made for its elasticity. If we 
compare these theoretical expressions of the 
nutational terms with the results of observations 
we shall notice at once some evident differences 
which cannot be ascribed solely to errors of 
observational data but may be explained as due 
to the dynamical effect of the earth’s core. 

We must keep in mind that equations of nuta- 
tion govern the motion of the earth’s angular 
momentum G. Since the position of an observa- 
tory is a definite place on the earth’s surface and 


TABLE I 


Nutation in longitude 


rived from 
Theoretical, observations, 
Term v sine We sine 
Main —6"869 sin Q —6"853 sinQ 
+0.008 cos 2 
Fortnightly —o"0812 sin 2¢ —0"0866 sin 2¢ 
} +0.0019 cos 2¢ 
Semiannual —o"507 sin 2© —0*533 sin 20 
Nutation in obliquity 
ivec from 
Theoretical, observations, 
=, @&— © 
Main +9%220 cos Q +9"198 cos Q 
9 —0.001 sin Q 
7 Fortnightly § +0%0884 sin 2¢ +0%0894 cos 2¢ 
; +0.0019 sin 2¢ 
~Semiannual +0%552 cos 20 +0"578 cos 20 
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an observer is thus always attached to the shell, 
so are the data obtained from astronomical ob- 
servations relevant of the motion of the shell 
alone but not of the earth as a whole. It follows 
that for a comparison with observations we 
should take equations governing the motion of 
the earth’s shell. We shall denote its angular 
momentum by G,. 

Since both G@ and G, are vectors lying in the 
equatorial plane XOY, they may be expressed 
by complex numbers, as follows: 


G = G (sine-¥ + is) (10) 
G, = G, (sin e-Wo + ié). (11) 


Let L be the couple arising from the attraction 
of the moon and sun on the earth’s shell. It is 
easy to show that 


L = hG, (sin e-} + 78), (12) 


where h is the ratio of H., the mechanical ellip- 
ticity of the shell to that of the earth as a whole. 
Since we deal now with the shell, the effect of 
the core should be considered as an action of ex- 
ternal forces. Denoting the moment of these 
forces by M, we may write an equation for the 
motion of G, in the following form: 


G,=L+ M. (13) 


The couple M transfers the angular momentum 
between the shell and the core but does not 
affect the momentum of the earth as a whole. 


Let us put 
M=X +1yY. (14) 


Substituting (11), (12) and (14) in (13), we shall 
have 


7. ae 
= G, [sin € (Wa a hy) a5 1 (€, ae he) }. (15) 


Let us compare now this expression for M with 
that obtained on the assumption that the core is 
rigid or, in general, that no motion of the core 
relative to the shell is possible. For this special 
case we denote the couple arising from a mutual 
action between the core and shell by 


M' = X'+1Y’, (16) 
and the angular momentum of the shell by G,’. 
Its direction is practically the same as that of 
the vector G and, consequently, its motion is 
governed by the equation 


G,’ = G, (sin e-¥ + 78), (17) 


in which y and « are the same as in (10). 
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Substituting the values of G,’ and L given by 
(17) and (12) in the following equation 


M’ =G,' — L 
we find 
X'+ iY’ =G, (1 —A)(sine-¥ + ze). (18) 
Let 
A, = the equatorial moment of inertia of the 


core, 


A = that of the earth as a whole, 
H, = the mechanical ellipticity of the core. 
According to K. Bullen (1936) 
A./A = 0.112, H, = 0.0026. 
Then 
h = 1.027. 


If we denote the earth’s angular velocity by n 
and put 


= wyt, 


2 (C = Wel, 20 ast 


we shall have 
@1 = — 0.000146 n, 


w2 = + 0.07300 n, 
w3 = + 0.00547 n. 


The substitution of theoretical expressions for y 
and ¢€ — ein (18) leads to the following equation: 


X’+iY’=U,4+ U,+0,+ WM (19) 
‘in which _ 
Uy =+ 072170:G,erimt, 
U2 =— 0%032u,G,e7 
(20) 


Vi =+ 070023w0G,et it, |” 
Wy =+ 0%014wsG,er iat 


NUTATION AND THE VARIATION OF LATITUDE 


By HAROLD JEFFREYS 
St. John’s College, Cambridge, England 


_ Abstract. A theoretical discussion by the author and R. O. Vicente uses geophysical estimates of the mechanical prope 
ties of the earth’s shell and two extreme models for the core, chosen to make the mass and moment of inertia corre 
period found for the variation of latitude is in good agreement with observation. The 18.6-year nutation is in b 
agreement than has been found previously but is still not altogether satisfactory. 


It is well known that the period of the 14- 
monthly variation of latitude is greatly affected 
by the elasticity of the earth. For a rigid earth 
the period would be about 305 days. The actual 
period is rather uncertain but can be taken as 
440 days with an extreme uncertainty of 15 
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Now let us take the values of Wa and eg — €o fro 
our Table. Being substituted in (15) they giv 


X + iY = (1.09 +/0.021) 0; ; 
+ (1.09 + 0.13?) U2 a 
+ oe 0.43 ae 0.832) Vi — 0.86W3. ( 


Having regard to the uncertainties of the 
served values as well as the computed values’ 
N and h, it is not easy to say to what extent th 
result is trustworthy. Nevertheless I should 
to point out the following conclusions which, 4 
my opinion, deserve some consideration : 


(1) The actual magnitude of the vector Us 
inferred from observational data, is larg 
than that obtained theoretically for 
rigid core; : 
The actual directions of the vectors | 
and W, are opposite to those for the ri 
core. 


(2) 


At first glance these conclusions seem to c 
tradict one another. However this contradictie 
vanishes under more close consideration. The ob 
served changes of the vectors mentioned abot 
agree in sign with that which would be expecte 


core. A quantitative comparison makes no se 
because of the lack of accuracy of observati 


data. a 
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days. The difference was for a long time our bes) 

datum on the elasticity of the earth as a wh 

Seismology has shown the earth to hav 

' central core, with a radius of about 0.55 of t 

‘of the outside; this does not transmit transverst 
waves and is presumably liquid. Seismology ha 
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lso led to rather detailed information about the 
,iy, lastic properties of the shell; if the rigidity at 
jae base of the shell was continued to the center 
jae elastic yielding in the variation of latitude 
rould be too small and the calculated period too 
nort, but if the core is liquid the yielding seemed 
kely to be enough. 

he _ The perturbations of the earth by the moon 
lei vive an estimate of the mass of the moon. Given 
iigine rate of precession and this mass the ratio 
i “€ — A)/C can be found, C and A being the prin- 
icy -ipal moments of inertia. There are many periodic 

“ouples due to the sun and moon, tending to 
; hift the earth’s axis, and their effects are called 
, \utations. The chief period is 18.6 years and 
‘d “ilts the axis with an amplitude of 9’207 + 
7 "003. The nutations can, however, all be cal- 

ulated from the mass of the moon and (C — A) 

_C, on the hypothesis that the earth is a rigid 
‘ody. The theoretical value is close to 9”2271. 
:. Jackson (1930) was the first to insist that the 
0 - ifference is genuine. 
jc | Poincaré (1910) studied the effects of a fluid 
te ore and elasticity of the shell. He showed that 
|) lasticity could have little effect on the 18.6-year 
pe ~ utation, but that the effect of a liquid core 
\y aight be considerable, as indeed Kelvin and 
)y thers had already found. He had, of course, no 
iy juantitative data, and his work was largely over- 
-ooked, though part of it is reproduced in Lamb’s 
_Tydrodynamics. Application to the earth 
Jeffreys 1948), assuming a rigid shell and a 
‘quid core of the actual dimensions, showed that 
he fluidity of the core would shorten the free 
ih yeriod from about 305 to 270 days, leaving more 

or elasticity to do to bring it up to the observed 
talue. The 18.6-year amplitude would be re- 

-luced by about 1 part in 150, which is in the 

ight direction but too great. I found, however 
Jeffreys 1949), that both effects would be 

_ educed by elasticity of the shell. 

The problem of the elastic yielding of the 
arth’s shell is not very difficult if it is treated as 
-iomogeneous and incompressible. With the 
_\ctual distribution of density and elasticity it is 

if yi very difficult, involving the numerical solution of 
/hree differential equations of the second order. 
it was solved by H. Takeuchi (1950) who used 
wo: models proposed by Bullen, which cannot be 
i tj lar wrong. 

R. O. Vicente and I undertook the allowance 
or the core in 1954; we made independent solu- 
ions but the results were in good agreement. 
Cakeuchi’s solution was a statical one, expressing 


baal! 


ah 
7 


L 
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all displacements of the shell in terms of six 
boundary values, but his relations will hold for 
periods of a day or longer. For the core it was 
necessary to retain the inertial terms. We used 
two models. In both, the mass of the core and its 
moments of inertia were made to agree with 
Bullard’s theory of the figure of the earth. In 
one (Jeffreys and Vicente 1957a) the core was 
treated as homogeneous and incompressible, with 
an extra particle at the center. In the other 
(Jeffreys and Vicente 1957b) it was taken as 
having a distribution of density of the type 
p = ko — kir*/ay’, the variation of density being 
taken as wholly due to compression. The truth 
should lie about midway between these models. 

The method was to introduce pairs of La- 
grangian coordinates to represent (1) a general 
rotation, displacing the axis, (2) normal displace- 
ments at the outer and inner boundaries, (3) a 
displacement of the core, linear in the coordi- 
nates, (4) and (5) two pairs of cubic terms in the 
coordinates for the core. These are of sufficient 
generality to represent any displacement of the 
core up to cubic terms in the coordinates, and all 
shell effects are represented by (1) and (2). The 
problem was thus reduced to a Lagrangian one of 
small oscillations, and the rest was arithmetic. 

For both models the free period was found to 
be 392 days, but allowance for the ocean would 
raise this to about 430 days. The correcting 
factor for the 18.6-year nutation was —0.0036 
for the central particle model and —o.oo11 for 
the Roche model; thus, compressibility of the 
core makes quite a large difference. The most 
probable value is about —0.0022, which is 
about right. Solutions were also made for the 
fortnightly and semiannual nutations, and the 
bodily tide numbers were calculated for the cor- 
responding tides. 

The theory predicts that the effect in longitude 
would be greater than that in obliquity. So far 
Fedorov is the only astronomer to have tried to 
estimate them separately from observation. In 
his announcement to the I.A.U. in 1953 he found 
the effect in obliquity to be the larger, but in his 
memoir just published he no longer maintains 
this. This is a most valuable piece of work. (A de- 
tailed comparison will appear in the Monthly 
Notices.) 

I very much doubt whether motions in the 
core can have any important effect on the 
general rotation of the earth. Bondi and Lyttle- 
ton found that there is much difficulty in explain- 
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ing how the couple due to tidal friction can be 
transmitted to the core by viscosity, and there 
would be more difficulty for periods from 0.5 
year to 3000 years. Electromagnetic coupling 
has been suggested, but so far as I know there 
has been no quantitative examination. 


LES RELATIONS ENTRE LES MOUVEMENTS DU POLE ET LES 
FLUCTUATIONS DE LA VITESSE DE ROTATION DE LA TERRE 


Par PAUL MELCHIOR 
Observatoire Royal de Belgique, Bruxelles, Belgique ee 
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Résumé. On expose comment on peut tenter une ample synthése des résultats d’observations astronomiques telles qu 
celles du mouvement du pdle et des fluctuations de la vitesse de rotation dela Terre avec des observations de caracteri 
purement géophysique: les marées terrestres, la circulation atmosphérique générale, le retrait progressif des glaciers, etc 

On ne peut espérer donner dés a présent une conclusion 4 chacun des problémes soulevés par la confrontation de ce 
données si diverses et l’on relévera plusieurs contradictions. L’effort accompli dans le cadre de l’Année Géophysique per 
mettra peut-étre de lever certaines d’entre-elles et de parfaire la synthése évoquée ici. 


Se proposer l’étude des anomalies que nous 
révéle l’observation de la rotation de la Terre, 
c’est aborder un probléme de géophysique avec 
des méthodes d’investigation astronomiques. 
Qu’il s’agisse de l’étude des composantes du 
mouvement du pdle, de l’étude des fluctuations 
de la vitesse de rotation du globe ou de prob- 
lémes plus particuliers comme celui de la con- 
stante de nutation, l’on devra faire appel a des in- 
formations purement géophysiques (séismologie, 
marées terrestres notamment) si l’on veut pro- 
gresser dans l’interprétation des phénoménes 
observés. Géophysiciens et astronomes étudient 
‘les mémes propriétés physiques de la méme 
planéte par des voies différentes. Une synthése 
riche en perspectives se prépare dans le cadre 
d’une doctrine commune que par analogie avec 
l’astronomie géodésique nous appellerions volon- 
tiers ‘‘astronomie géophysique.”’ 

Nous diviserons notre exposé sur les relations 
entre les mouvements du pdle et les fluctuations 
de la vitesse de rotation de la Terte en trois 
parties qui sont trois problémes essentiellement 
différents au point de vue du géophysicien mais 
qui, dans l’observation astronomique, sont 
étroitement associés: ; 


Premier probléme. Relation entre la composante 
libre du mouvement du pdle, les fluctuations a 
période semi-mensuelle, mensuelle et semi-an- 
nuelle de la vitesse de rotation de la Terre, les 
marées terrestres, la séismologie, le relévement 
isostatique. 

Deuxiéme probleme. Relation entre la com- 
posante forcée annuelle du mouvement du péle, 
la fluctuation a période annuelle de la vitesse de 


rotation de la Terre, la météorologie et l’océano. 
graphie. 

Troisiéme probléme. Relation entre le déplace. 
ment séculaire du péle, les variations séculaires 
de la vitesse de rotation de la Terre, la glact- 
ologie. q 

I 


1. L’étude des phénoménes de marées ter. 
restres et l’étude du mouvement libre du pdk 
instantané de rotation a la surface de la Terre 
rentrent dans le cadre d’une seule et méme 
théorie. 

En effet, si la Terre était une Terre idéale 
indéformable comme le suppose la Mécanique 
Céleste classique, on pourrait rendre exactement 
compte par le calcul des caractéristiques de ces 
deux phénoménes et méme en établir des éphé- 
mérides. Mais les forces qui entrent en jev 
produisent dans les deux cas des déformations 
et engendrent des tensions dans tout le globe. 
Ces déformations a leur tour additionnant leur 
effet a celui du potentiel initial modifient les 
caractéristiques des phénoménes. En premiére 
approximation, comme il s’agit de phénoménes 4 
longue période vis a vis de la période de vibration 
propre de la Terre (1534™), on peut considérer 
que ces effets obéissent a la théorie statique. Ce 
sont donc les amplitudes des phénoménes qui 
sont modifiées: réduites ou augmentées suivant 
les manifestations du phénoméne que I’on con- 
sidére. 

Ainsi dans le cas du mouvement du péle, le 
pole d’inertie vrai tend a se rapprocher du pdle 
instantané de rotation et le suit dans son mouve- 
ment; comme toute observation astronomique 


— 


lu pole d’inertie vrai est impossible, nous sommes 
amenés a prendre pour origine de référence un 
vertain pdle d’inertie “moyen” que nous dé- 
‘inissons un peu arbitrairement (cf. Troisiéme 
Probléme) comme le barycentre de la polhodie 
‘};ouvrant six années (commensurabilité approxi- 
‘native des deux périodes fondamentales recon- 
-aues dans le mouvement du pdle). 

| Cette maniére inéluctable de procéder donne 
'apparence d’un allongement de période—c’est 
‘’allongement chandlérien—il n’est que l’image 
‘du déplacement du péle d’inertie vers le péle de 
rotation, déplacement dfi a I’action de la force 
centrifuge sur la matiére élastique qui constitue 
la Terre et dont l’amplitude est fonction du 
module de rigidité de la planéte. Tel est le lien 
existant entre la période de Chandler et les 
“propriétés élastiques de la Terre. 

| On voit qu’il s’agit en somme d’introduire dans 
“les calculs de la Mécanique Céleste les para- 
“métres représentant |’état physique interne de 
“la Terre. On avait dans le passé déja tenté de 
‘discuter l’effet de la distribution radiale des 
‘““‘densités (paramétre de Radau—aplatissement 
dynamique); On tente aujourd’hui d’introduire 
‘la distribution radiale du module de rigidité, le 
coefficient de viscosité, l’effet des convections 


— 


_ 2. Examinons d’abord de plus prés comment la 
“théorie permet d’associer 4 ce phénoméne astro- 
~nomique le phénoméne géophysique des marées 
_ terrestres. L’idée fondamentale de la théorie, due 
a Love, est que, le potentiel étant représenté 
_(avec une approximation suffisante) par une 
‘fonction harmonique sphérique du second ordre, 
toutes les perturbations produites dans la Terre 
‘par ce potentiel peuvent étre représentées par la 
‘méme fonction harmonique, affectée d’un facteur 
numérique d’amplitude correspondant a |’aspect 
du phénoméne particulier envisagé; ce facteur 
est appelé nombre de Love; pour certaines mani- 
festations de la déformation ce sont des combi- 
naisons algébriques simples des nombres de Love 
_ que l’on sera amené 4a introduire. 

_ A titre d’exemple, le déplacement radial U, la 
‘dilatation cubique D produits en tout point de la 
Terre par les forces dérivant d’un potentiel du 
‘second ordre W2 sont donnés par 


U = H(r) W2/g 
D = f(r) W2/g, 


r étant la distance du centre de la Terre au point 
considéré. Par voie de conséquence le potentiel 
dfi a la déformation et a la variation de densité 
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provenant de la dilatation cubique sera donné 
également par une expression de la forme 


V = K(r) We 
tandis que les déformations dans le méridien et 
le premier vertical s’écriront 
bW, 
g 60 


dW. 
g sin 0 6X 


L(r) L(r) 
Le potentiel W.2 peut étre aussi bien le potentiel 
d’attraction luni-solaire que le potentiel de la 
force centrifuge provenant du mouvement de 
l’axe de rotation par rapport a la Terre: 


We. = V — Vo = } wr? (sin? ¢’ — sin? d) 
ou ¢’ = ¢ — xcosA — ysini, 


(x y) coordonnées du Service International des 
Latitudes, \ = longitude. 


A la surface de la Terre (ry = a) on pose 
H(@)=h IK @) = iE (ayia) 


qui sont les nombres de Love usuels: 


h = le rapport de la hauteur de la marée ter- 
restre a la hauteur de la marée statique 
océanique correspondante, 

k = le rapport du potentiel additionnel dfi a 
cette déformation au potentiel déformateur 
lui-méme. 


L’amplitude des déviations de la verticale par 
rapport a la crofite, dues aux marées, est affectée 
du coefficient y = 1 + k — h, les déviations de 
la verticale par rapport a l’axe du monde du 
coefficient A = 1 + k — 1, les variations de |’in- 
tensité de la pesanteur, du coefficient 6 = 1 +h 
— 32k. On a trouvé y = 0,706; 6 = 1,20; 
A = 1,15, d’ot h = 0,482 et k = 0,188. 

Pour fixer les idées sur les valeurs numériques 
de ces nombres, notons que: 


1° pour une Terre h=1 k = 0,6 10:3 
fluide homogéne 

2° pour une Terre in- h=0 k=0 Ll=0 
finiment rigide 

3° pour la Terre réelle = 0,482 k =0,188 1 = 0,06 


Les nombres h, k, / sont évidemment liés a la 
distribution radiale des constantes élastiques et 
de la densité du milieu par des équations diffé- 
rentielles (Herglotz) complexes (généralisation de 
’équation de Clairaut) dont l’intégration nu- 
mérique a été exécutée pour la premiére fois par 
Takeuchi et ensuite par Molodensky. Toute une 
série de modéles de la constitution interne du 
globe ont ainsi été éprouvés et certains donnent 
des résultats assez voisins de l’expérience bien 


88 


qu’un accord parfait entre les nombres y et 6 ne 
soit pas encore obtenu. L’introduction de la 
graine dans les calculs le permettra sans doute. 

On montrera maintenant comment la valeur 
numérique du nombre k& détermine la valeur de 
la période du mouvement libre du pdle instan- 
tané de rotation a la surface de la Terre, dite 
période de Chandler. 

Le potentiel perturbateur W. di au mouve- 
ment du pdle (force centrifuge) peut étre re- 
présenté par une fonction harmonique tesserale 
du second ordre, tout comme le potentiel des 
marées diurnes (les marées semi-diurnes et a 
longue période étant représentées par des fonc- 
tions harmoniques du second ordre respective- 
ment sectorielles et zonales). 

La déformation consécutive a l’action de ce 
potentiel W»2 et qui lui est proportionnelle induit 
un nouveau potentiel additionnel kW». et intro- 
duit en méme temps de petits produits d’inertie 
non nuls dans les équations d’Euler. Ces produits 
d’inertie contiennent donc en facteur le nombre 
k et peuvent aisément €tre exprimés en fonction 
de la différence des moments principaux d’inertie 
polaire et équatorial C — A: 


=—kyu(C— A) 


— kxp(C— A) 


D 
joe 
Ff =0 


avec 


oti (x y) sont les coordonnées du Service Interna- 
tional des Latitudes, 


vitesse angulaire de rotation de la Terre, 
= rayon de la Terre, 

intensité de la pesanteur, 

aplatissement de la Terre. 


mm RE 
| 


Les équations du mouvement s’écrivent alors 


dx 6 — A 
aot oe A Fe) oy 0, 
d C-A 
ie a (t — ku) wx =0, 


et la période du mouvement libre est 


relation qui est indépendante de toute hypothése 
concernant la structure interne de la Terre, sauf 
dans la définition et la valeur de b. 

On voit aussit6t que pour k = 0, cas de la 
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Terre rigide, on retrouve la période classiqu 
d’Euler: 


20 A 
wo C—-A 


b= (1 -%). 
Ve Ti 


La valeur trés sfire donnée pour & par la discus 
sion des observations de marées terrestre 
(k = 0,188) conduirait 4 une période de Chand 
ler de 382 jours, en net désaccord avec tous le 
_résultats d’analyse harmonique effectués sur le 
courbes décrites depuis 70 ans par le péle. 

Notons toutefois que l’on devrait apporte 
une correction 4 cette période pour tenir compt 
de la contribution des océans dans les produit 
d’inertie de déformation. On arriverait ainsi ; 
une valeur maximum de k d’environ 0,23 soit 
une période de 400 jours: le désaccord ave 
l’observation est a peine atténué. 

3. L’analyse harmonique de toutes les donnée 
récoltées par le Service International des Lat 
tudes depuis sa fondation fournit en effet | 
valeur moyenne de 429,5 jours pour la période d 
Chandler (Melchior). Cette analyse harmoniqu 
a pu étre menée de bien des maniéres mais le 
résultats ne sont pas significativement différents 

Cependant, si l’on admet A priori que | 
mouvement obéit 4 une certaine loi analytiqu 
et si l’on tient compte de la forme de cette Ic 
pour exécuter le calcul, on peut étre conduit | 
des résultats qui ne sont pas nécessairement € 
accord avec ceux des analyses exécutées san 
schéma précongu. C’est ce qu’ont fait Walker e 
Young en admettant, aprés Dyson et Jeffrey 
que le mouvement est amorti dans les condition 
établies par le postulat de 1|’élastico-viscosité d 
Maxwell. Dans ce cas les équations s’écriven 
sous la forme générale: 


— &; + Axi + xx; = 0. 


En appliquant le schéma d’autorégression d 
Yule, Walker et Young parviennent a une périod 
de 440 jours (désaccord encore plus prononcé ave 
les marées terrestres) et un temps de relaxatio: 
2m 
oe 
_ plitude @ est réduite 4 ea soit 0,37 de l’ampli 
tude initiale. 

Gutenberg a mis en doute ce résultat, mon 
trant que si l’on admet pour le module de rigidit 
les valeurs déduites des vitesses de propagatio: 


= 305 jours sidéraux. 


On a donc 


de 3 ans. r = la durée aprés laquelle l’am 


1959 April 


des ondes sismiques confrontées avec l’amplitude 
des marées terrestres et en admettant pour le 
coefficient de viscosité la valeur déduite de la 
vitesse du relévement finno-scandien (Niskanen), 
on arriverait 4 un temps de relaxation de plus- 
‘ieurs centaines d’années. 

Aucune explication n’a d’autre part pu étre 
trouvée pour un amortissement aussi fort. 
Melchior a mis en doute sa signification en 
montrant que les formules de Yule reposaient 
‘sur le postulat de Maxwell, lequel n’a aucune 
base physique: on peut donc douter de la vali- 
dité de ces équations adoptées pour représenter 
le mouvement du pdle et penser que les valeurs 
numériques trouvées par Walker et Young sont 
des paramétres donnant une bonne représenta- 
tion mathématique mais dénués de sens physique. 
__ Cependant une difficulté supplémentaire réside 
‘dans le fait que la période du mouvement du 
‘pole apparait comme nettement variable et 
qu’aucune démonstration théorique n’a encore 
pu étre donnée pour expliquer un tel phénoméne. 
_ 4. Un inconvénient majeur réside dans le fait 
“que la suite des coordonnées (x y) du péle man- 
que d’homogénéité pour des raisons diverses que 
‘nous énumérons briévement: 
| 


a i 


ee ee 


° La nécessité d’assurer dans chaque 
ine d’étoiles observé une compensation aussi 
parfaite que possible des mesures micrométriques 
(de sorte que les erreurs inévitables sur la valeur 
admise pour les pas de vis ne se répercutent pas 
sur les latitudes déduites des observations) 
oblige a des révisions périodiques du programme 
d’étoiles A observer: substitutions de certaines 
-étoiles ou paires d’étoiles rétablissant un équilibre 
que la précession tend a détruire systématique- 
ment au cours du temps. Ceci fut fait en 1906, 
1912, 1922, 1935, 1954. Mais l’effet de ces change- 
“ments de catalogue fut plus sensible encore 
lorsque les systémes auxquels ils étaient référés 
ont également été changés: AGKr en 1900, 
G. C. Boss depuis 1935. 
' 2° La méthode de réduction a varié selon les 
conceptions des savants qui furent chargés de la 
réduction d’ensemble des données. Les variantes 
“successives ont évidemment toutes été intro- 
‘duites dans un but de perfectionnement. Mais 
leur juxtaposition fait que le procédé actuelle- 
“ment: en usage différe sensiblement de celui qui 
i} fut mis en oeuvre au début du siécle par Albrecht. 
‘Les anciennes séries d’observations devraient 
donc étre revues selon les critéres actuellement 
admis. 


A titre d’exemples, on introduit maintenant 
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de petits termes de nutation (Ross) les termes 
d’aberration de Battermann, tandis que les com- 
posantes périodiques locales empiriques de 
Kimura ont été abandonnées. Par contre |’inter- 
prétation du terme z de Kimura et de |’erreur de 
fermeture est a présent bien assise. On peut 
logiquement considérer .le terme de Kimura 
comme dfi essentiellement aux erreurs de décli- 
naisons du type Aé, avec en surplus de petites 
contributions dues aux erreurs commises sur les 
valeurs admises pour les constantes fondamen- 
tales de nutation et d’aberration. 

L’erreur de fermeture est pour 60 pour-cent 
due a l’erreur sur la constante d’aberration 
(laquelle devrait €tre portée 4 20"51); pour le 
reste elle parait bien liée aux variations de la 
pression atmosphérique au cours de la nuit 
(Hattori). Il est bien évident que des effets 
météorologiques locaux assez difficiles a définir 
jouent un réle non négligeable dans ces deux 
caractéristiques des observations tandis que la 
nature de |’instrument, son type de construction 
ne joue pas de réle apparent (sauf pour le 
télescope flottant de Cookson, qui est assez 
sensible aux effets du vent, ce qui n’est guére 
surprenant). 

Ainsi donc une nouvelle réduction d’ensemble 
doit avant tout se baser sur un catalogue homo- 
gene établi dans le systéme du FK4 (la con- 
struction de ce catalogue est en cours a 1|’Obser- 
vatoire Royal de Belgique). Elle doit étre faite 
sur la base d’une méthode de réduction soigneuse- 
ment mise au point, méthdde qui ne différerait 
pas sensiblement de celle du Professeur Cecchini, 
sauf considération aux variations a longue 
période des latitudes locales qui devraient étre 
éliminées selon les vues d’Orlov mais avec une 
combinaison d’ordonnées plus sélective que la 
sienne (cf. III.1). 

Dans cette réduction, on devra adopter la 
valeur 20"51 pour la constante d’aberration et 
une valeur améliorée de la constante de nutation. 
Enfin, l’effet des marées terrestres sur la direction 
de la verticale locale par rapport a l’axe du monde 
doit étre pris en considération puisque l’opération 
inverse a permis de déduire la marée des observa- 
tions méme de latitude du Service International 
et d’en donner une définition des nombres de 
Love qui est trés cohérente avec les observations 
purement géophysiques (cf. I.6). 

Nous pensons que ce n’est qu’aprés un tel 
travail—a vrai dire trés long et trés délicat—que 
l’on pourra avoir la certitude que les caractéris- 
tiques des composantes du mouvement du pdle 
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sont absolument indépendantes et du procédé 
d’observation et du procédé de réduction. 

5. Considérons maintenant les fluctuations de 
la vitesse de rotation de la Terre autour de son 
axe instantané de rotation. Parmi ces fluctuations 
il existe des composantes qui résultent des dé- 
formations de la Terre provoquées par des ma- 
rées terrestres du type zonal (et qui sont les 
marées dites a longue période). Les déforma- 
tions zonales modifient en effet le plus grand 
moment d’inertie C du globe et l’on peut s’at- 
tendre par conséquent a une variation périodique 
correspondante de la vitesse de rotation de la 
Terre en vertu de la formule: 


Ainsi l’on retrouvera l’effet de la marée solaire 
semi-annuelle Ssa et des marées lunaires men- 
suelle Mm (27,56 jours) et semi-mensuelle Mf 
(13,66 jours) dans les fluctuations de la vitesse 
de rotation. Elles ont été mises en évidence pour 
la premiére fois par W. Markowitz dans les 
observations au P.Z.T. L’amplitude des marées 
lunaires est de l’ordre de la milliseconde et 
contient en facteur le nombre k. En effet, le 
potentiel kW. engendré par la déformation 
zonale considérée s’exprime en fonction des 
variations du moment d’inertie principal C. 

L’amplitude de la composante périodique de 
-la variation de vitesse est donc bien liée au 
nombre k et par conséquent |’observation astro- 
nomique de haute précision est appelée a fournir 
bient6t une nouvelle valeur expérimentale de ce 
nombre. 

Cette méthode nous est trés précieuse en 
raison de l’impasse dans laquelle on se trouve 
actuellement par suite du désaccord entre la 
valeur ‘‘chandlérienne”’ de & et la valeur déduite 
des observations purement géophysiques de 
marées terrestres. 

En ce qui concerne la marée solaire semi-an- 
nuelle Ssa, N. Stoyko avait déjé montré en 
1950 qu’elle permettait d’expliquer la plus grande 
partie de la variation semi-annuelle de la vitesse 
de rotation. Les résultats de Markowitz, indi- 
quant une amplitude de 10 millisecondes, le 
confirment. 

6. Enfin, les méthodes astronomiques mises en 
oeuvre pour l’étude de la rotation de la Terre 
constituent encore deux méthodes indépendantes 
pour la détermination expérimentale des nom- 
bres de Love. On sait, en effet, que les déviations 
de la verticale sont observées par les géophysi- 


~mouvements de masses d’air et d’eau obéissam| 


‘ 
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ciens a l’aide de pendules horizontaux tré 
sensibles placés a grande profondeur (de 1’ordre 
de 100 métres). Cette maniére d’observer le 
déviations de la verticale n’offre d’autre référ. 
ence que la crofite mobile elle-méme et introduit 
en facteur la combinaison y = 1 +k —h de 
nombres de Love. II est évident que ces dévia 
tions périodiques de la verticale affectent le: 
déterminations de latitude de haute précision e1 
que de trés longues séries d’observations telle: 
que celles du Service International des Latitude: 
(ou de certains observatoires s’étant consacrés é 
ce probléme) sont susceptibles de les mettre er 
évidence. Dans ce cas cependant les observation: 
ne se référent plus a la crofite mais a l’axe di 
monde et le coefficient d’amplitude est la combi 
naison A =1-+k-—J. La discussion de tou 
l'ensemble de données actuellement disponibl 
a permis de fixer pour ce coefficient la valeur 


A=1+k-—J=1,15 


qui est en excellent accord avec les mesure} 
géophysiques. 

Le second procédé astronomique de détermina 
tion expérimentale des nombres de Love résult; 
de la comparaison d’horloges a balancier di 
haute précision avec des horloges a quartz. Le 
horloges 4 balancier sont abandonnés pour I) 
conservation de l’heure mais leur marche pré 
sentait cependant cet intérét d’étre liée a Ij 
valeur instantanée de g, laquelle varie en fonctio] 
de l’attraction luni-solaire et de la déformatio] 
du globe qui en résulte; les horloges a quartz ej 
sont indépendantes. Le facteur d’amplitud| 
rendant compte de l’effet de la déformatio; 
élastique est ici, comme pour les gravimétres, |) 
combinaison 6 =1+h4-— 4k. Une étude tré 
soigneuse due a N. Stoyko et portant sur quat 
années et cinq pendules de 1l’Observatoire dj 
Paris a donné pour résultat 


6 = 1,21 3 


en accord parfait avec les résultats oe 
triques: 6 = 1,20. / i 
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1. L’existence d’une composante annuell| 
importante aussi bien dans les fluctuations de ]} 
vitesse de rotation que dans le mouvement d/ 
pole est évidemment liée aux déplacements d 
masse a période annuelle a la surface de la Terre} 


au cycle saisonnier (c’est donc un effet solair 
indirect). S 
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2. L’importance de cet effet dans le mouve- 
ment du pdle s’explique par l’existence d’une 
résonance avec la période propre (période de 
Chandler). 

Les équations du mouvement du pdle s’écrivent 


d ; : : : ; 
Ae tin) — ix @ +49) =— x (E+in), 


ou (x y) sont les coordonnées du pdle instantané 
de rotation publiées par le S.I.L., (& 7) les co- 
ordonnées du péle instantané d’inertie, 7 = 27/y, 
période de Chandler, T = 27/a, période annuelle. 
On les intégre aisément en admettant le mouve- 
ment du pdle d’inertie elliptique (le calcul 
daprés les données météorologiques montre 
qu'il est quasi linéaire) : 


E + in = next + n'e~iat 
et l'on a 
eter eke? — mm. e864" 
Re'x' représente le mouvement libre chandlérien 
(circulaire) tandis que les axes du mouvement 


elliptique annuel du péle instantané de rotation 
sont donnés par 


b ’ n nN 
=m —m= _ 
r+ T— 1 
ae n' n 
a=m m= 
r+! 7—I1 


ie en (a, 7) et (6,7) qui sont trés bien 
satisfaites par les résultats de l’analyse harmoni- 
que des mouvements du pdle (1900-1955) et qui, 
compte tenu des variations observées de la 
période 7, conduisent pour les axes de l’ellipse 
décrite par le péle d’inertie aux valeurs: 


A=n-+7n'=o0'012, B=n-—n' = 0’040 

2n trés bon accord avec les calculs géophysiques 
de Rosenhead (calcul de l’effet des déplacements 
je masses d’air et d'eau, compte tenu d’une com- 
densation due a la flexion de la crofite terrestre) : 


A=olo1r B= 0%044. 


3. Ce raisonnement suppose que l’amortisse- 
nent du mouvement libre n’est pas sensible a 
‘échelle de la précision actuelle, c’est a dire que 
e temps de relaxation serait supérieur a 20 ans. 
Mais en admettant avec Walker et Young un 
‘emps de relaxation de trois ans, la résonance 
igissant sur le mouvement forcé annuel est 
ortement réduite et si l’on prend aussi pour la 
%ériode de Chandler la valeur de 440 jours, 
‘amplitude du mouvement du pdle d’inertie 
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devrait étre multipliée par un facteur trois pour 
rendre compte de la composante annuelle ob- 
servée dans le mouvement du péle de rotation. 
Les axes de l’ellipse représentant le mouve- 
ment du pdle d’inertie devraient donc avoir pour 
longueur 
Ag 07036" 5 —=.07,120 


ce qui ne s’accorde plus avec les calculs de 
Rosenhead. 

Cependant Munk et Groves ont pris en con- 
sidération l’effet des vents qui se manifeste par 
un échange de quantités de mouvement et une 
pression sur les chaines de montagnes (l’Him- 
alaya, orienté selon un paralléle, jouant 14 un 
role essentiel) ; ils ont obtenu: 


At —tOAOlie > — OalTA 


mais l’incertitude qui régne sur leurs données de 
base et précisément sur l’effet dans |’'Himalaya 
qu’ils ont traité d’une maniére hypothétique 
rend ces valeurs peu sfires. 

Par contre si l’on admet que l’amortissement 
n’est pas sensible (r > 20 ans), et que les vari- 
ations observées dans la période de Chandler sont 
réelles (bien qu’inexpliquées), il devient aisé 
d’expliquer les variations d’amplitude corres- 
pondantes de la composante annuelle comme 
une conséquence de la variation de l’effet de 
résonance. 

Ceci montre en tous cas que le probléme n’est 
pas résolu et que les résultats météorologiques 
de 1’A.G.I. permettront certainement de le 
traiter avec plus de sécurité. 

Notons encore que si le mouvement subissait 
un tel amortissement, il faudrait encore trouver 
la source de son entretien; on n’a pu la trouver 
ni dans les phénoménes atmosphériques irrégu- 
liers ni dans le noyau liquide. 

Quoiqu’il en soit notre conclusion est que la 
précision actuelle (1900-1958) des détermina- 
tions de la position du péle ne permet pas une 
détermination du temps de relaxation si celui-ci 
est supérieur a 20 ans, ce qui parait bien étre le 
cas. 

4. Les considérations qui précédent évoquent 
immédiatement |’existence d’une corrélation avec 
la fluctuation annuelle de la vitesse de rotation 
de le Terre. Cependant dans ce cas les déplace- 
ments de masses ne jouent qu’un faible réle et 
ce sont surtout les vents qui, par des échanges de 
quantités de mouvement, expliquent le phé- 
nomeéne. 

Le léger retard par rapport a la phase des 
saisons astronomiques n'est pas inattendu, 
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TABLEAU I. LES CORRELATIONS AVEC LA FLUCTATION ANNUELLE 
(PHASE t EXPRIMEE EN FRACTIONS D’ANNEES) 


Vitesse de rotation* minimum moyenne maximum moyenne Phase 
Accélération nulle maximum(-+) nulle maximum(—) moyenne 
0,104 0,354 0,604 ‘0,854 0,4790 
Elongation du pdle d’inertiet minimum maximum minimum maximum 
Elongation du pdle de rotation maximum minimum maximum minimum 
0,102 0,357 0,612 0,854 0,4812 
Milieu des saisons hiver printemps été automne 
astronomiques 0,096 0,342 0,597 0,841 0,4702 


*D’aprés N. Stoyko. 
t D’aprés P. Melchior. 


puisque les phénoménes considérés sont dus a 
une accumulation ou a une dissipation d’énergie 
par l’atmosphére et les océans. 


III 


1. A ces phénoménes annuels, chandlériens et 
luni-solaires, se superpose aussi bien dans le cas 
du mouvement du pdle que dans le cas des fluctu- 
ations de la vitesse de rotation de la Terre un 
effet dit séculaire en ce sens qu’on n’y reléve 
pas de périodicité définie soit parce que nos ob- 
servations ne s’étendent pas suffisamment loin 
dans les siécles passés, soit parce qu’il n’en existe 
pas. 

Ces dérives ne sont pas linéaires, elles présent- 
ent des inversions, des perturbations, des dis- 
continuités. Comme dans les deux cas précédents 
on a tenté naturellement d’établir une corrélation 
entre les accidents qui ont affecté la vitesse de 
rotation et ceux qui ont affecté la position du 
pole. 

2. La question d’une dérive séculaire possible 
du pdle est d’une trés grande importance en 
géologie et en géophysique, mais pour |’instant, 
nous ne pouvons nous baser que sur a peine plus 
d’un demi-siécle d’observations réalisées par le 
Service International des Latitudes. Les observa- 
tions antérieures ne sont ni de la méme précision 
ni de la méme homogénéité; en outre elles sont 
basées sur des catalogues d’étoiles.de précision 
nettement insuffisante. Méme dans le cas du 
Service des Latitudes, l’homogénéité est loin 
d’étre parfaite et les discontinuités dont nous 
avons parlé plus haut (I.4) empéchent, comme 
nous le montrerons plus loin, de tirer actuelle- 
ment une conclusion tant soit peu sfire sur cette 
dérive hypothétique. Le probléme du mouve- 


ment séculaire est étroitement lié A la définition — 


exacte de la position du pdle d’inertie moyen, de 
la latitude moyenne de chaque station et de la 
stabilité de cette latitude. 
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3. Le probléme de la détermination de la lati 
tude moyenne des stations du Service des Lati 
tudes est un probléme controversé du point d 
vue du procédé de réduction a adopter pour | 
observations de latitude. II sera nécessaire d 
trancher la question au moment ot l'on s 
proposera d’exécuter la révision d’ensemble d 
toutes les observations faites depuis 1900. On n 
peut, en effet, prétendre que les stations de lati 
tude soient douées d’une stabilité parfaite 
Certaines d’entre elles présentent méme de 
glissements relativement importants (Mizusawa 
Tschardjui). Toutes les stations, comme 1’: 
montré A. J. Orlov, présentent en outre de lente 
fluctuations dans leur latitude moyenne, fluctu 
ations dont l’origine n’a pas été définie mai 
d’oti un mouvement réel du péle doit étre exclu 
car on constate qu’elles sont différentes pour de 
stations-situées sur le méme méridien ou mém 
pour deux instruments opérant en une mém 
station. Il faudra donc éliminer ces fluctuation 
de la latitude moyenne du calcul de la polhodie 
sous réserve de les discuter séparément. 

Or, tous les calculs du Service des Latitude 
sont effectués au milliéme de seconde, la latitud 
moyenne devra donc étre donnée pour chaqu 
dixiéme année avec cette précision par un 
combinaison adéquate des données encadran 
ce dixiéme d’année (méthode de Labrouste, 
Orlov a proposé que l’on élimine ces effet 
locaux des variations de latitude avant le calet 
des coordonnées du pdle en définissant 4 tou 
instant la latitude moyenne par la combinaiso 
linéaire d’ordonnées, 


4 
om = = dX (b: + ders + dite + e411). 
t=0 


_ Il est évident qu’on ne pourra admettre pour le 
_réductions des latitudes moyennes présentar 


d’une dixiéme d’année a l’autre des variations 
caractére brusque ou accidentel (méme dans ] 
milliéme de seconde), c’est 4 dire que la courb 


jeprésentant les fluctuations de la latitued 
moyenne dont on tiendra compte doit étre trés 
idoucie. 

L’application de la méthode d’Orlov montre 
jue cette combinaison d’ordonnées ne filtre pas 
‘igoureusement les courtes périodes. Or. tel 
levait étre l’objet de l’opération. La combinaison 
TOrlov a été construite, en effet, de maniére a 
yrésenter, dans sa courbe de sélectivité un zéro 
our m = 10, période annuelle, et un zéro pour 
1 = 12, période de Chandler, (” étant compté en 
lixiémes d’années). Or ceci n’assure pas |’élimi- 
vation rigoureuse des effets dfis au mouvement du 
6le car les périodes des composantes de celui-ci 
ie sont pas définies avec cette sfireté, en particu- 
lier la période libre. 

On est en présence, a cet égard, de deux 
typothéses: ou bien l’on a affaire 4 une période 
ibre oscillant entre les valeurs ” = 11,3 et 
+ = 12,0 ou bien I’on a affaire A un spectre de 
andes correspondant a la résonance du systéme 
“nécanique excité par des fluctuations désordon- 
‘1ées (d’origine non précisée) ce spectre de bandes 
‘st assez large. Pour éliminer ces composantes, 
‘| importe peu de fixer exactement les zéros de la 
jourbe de sélectivité 4 m = 10 et m = 12; il est 
»lut6t nécessaire de choisir une combinaison 
‘omplexe d’addition qui posséde une série de 
‘éros bien répartis entre » = 3 et m = 15, de 
‘naniére a filtrer rigoureusement toutes les 
ourtes périodes. Du point de vue purement 
‘héorique on peut donc critiquer la méthode 
‘Orlov parce qu’elle ne répond pas exactement 
Ux caractéristiques physiques du mouvement 
‘lu pdle. En outre, sa sélectivité n’est pas suff- 
‘ante pour assurer |’élimination des courtes 
»ériodes (définies par ” < 15). 

_ Pour choisir une méthode plus sélective, nous 
‘ous sommes donné pour conditions de centrer la 
-ombinaison sur un dixiéme d’anneé exact et 
‘’éliminer toutes les périodes inférieures An = 
5. La combinaison complexe de rang supérieur 
» introduire sera telle que n = 15 >2m+1 
‘oit m <7; cette combinaison sera s;. Nous 
ons porté notre choix sur la combinaison 
4855657 (notations de Labrouste). 

Le facteur d’amplification pour les longues 
€riodes est 20 dans la combinaison d’Orlov et 
lans la combinaison $455s657, 99 X 195 = 19,305. 
es courbes de sélectivité montrent évidemment 


En particulier la faible oscillation de période 
= 0,7 an, qu introduit la combinaison d’Orlov 
ans les variations de latitude et qui correspond 
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au maximum que présente la courbe de sélecti- 
vité pour cette période, est évidemment inexis- 
tante avec la combinaison s4555¢57. 

4. Ceci étant posé, considérons la maniére 
dont le mouvement séculaire du pdéle ‘‘moyen”’ a 
pu étre mis en évidence. Sekiguchi a opéré des 
moyennes glissantes sur six années sur les co- 
ordonnées du Service des Latitudes x, y; il a da 
cependant corriger a priori les coordonnées 
publiées des discontinuités qui ont été introdui- 
tes dans le systéme de référence a chaque change- 
ment de catalogue. 

A la date d’un changement de catalogue (cf. 
1.4), se produisent en effet a la fois une discon- 
tinuité et un changement de direction et de 
vitesse de la dérive. La discontinuité se traduit 
par une translation du systéme d’axes de co- 
ordonnées qui résulte en fait d’une variation 
inégale des latitudes moyennes des diverses 
stations; ces latitudes moyennes sont calculées 
par le Bureau Central sur chaque intervalle 
correspondant a un méme catalogue. 

Sekiguchi a utilisé des transformations de 
coordonnées calculées par Kimura et Nicolini 
d’aprés les changements apportés aux latitudes 
moyennes adoptées pour chaque station. La 
validité de ces transformations a été critiquée 
avec de bonnes raisons par Cecchini (déplace- 
ment réel de stations, telle Mizusawa). Pour 
déterminer sans hypothése préalable et les 
dérives et les translations du systéme d’axes 
Melchior a appliqué aux coordonnées originales 
(x y), sams aucune correction préalable et séparé- 
ment sur chaque intervalle homogéne, la com- 
binaison d’ordonnées 4555657 dont les propriétés 
ont été décrites plus haut. Les dérives trouvées 
sont généralement linéaires et l’on peut se 
permettre de les extrapoler jusqu’aux frontiéres 
de chaque intervalle considéré. La discontinuité 
aux frontiéres donne la transformation de co- 
ordonnées au moment du changement de pro- 
gramme par une voie tout a fait indépendante 
de la précédente et l’on constate que les formules 
de Kimura-Nicolini sont valables. Ceci implique 
que les stations, telle Mizusawa ne se déplacent 
pas réellement et qu’il s’agit de dérives au sens 
d’Orlov, c’est a dire fictives et dépourvues de 
sens astronomique ou géophysique. 

Ces dérives en x et en y sont absolument différ- 
entes d’un programme a l'autre. En outre Phillip- 
pov a démontré que les longues séries d’observa- 
tions de Greenwich, Poulkovo et Washington 
sont mieux représentées par une polhodie sans 
dérive que par la polhodie avec dérive. Nous 
pensons donc qu’on ne peut rien conclure pour 
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l’instant au sujet de la dérive du pdle moyen. 
Comme l’a remarqué avec pertinence Fedorov, 
si nous sommes forcés d’abandonner pour |’in- 
stant toute tentative d’étudier le mouvement 
séculaire du péle, c’est évidemment moins grave 
que de parvenir a des conclusions fausses sur ce 
mouvement. 

Ainsi donc le mouvement séculaire devrait 
étre dissocié de celui de la polhodie, ce qui pré- 
senterait l’avantage de la centrer sur l’origine et 
de fournir des corrections de variations de 
longitude plus pratiques pour les observations 
horaires. Cela revient a éliminer les fluctua- 
tions de la latitude moyenne des diverses sta- 
tions (comme nous en avons déja montré la 
nécessité en IIJ.3) et les arguments qui s’opposent 
a l’introduction de stations indépendantes du 
S.I.L. dans le calcul de la polhodie elle-méme 
(et qui reposent essentiellement sur la nature du 
terme annuel z de Kimura) ne jouent pas dans le 
cas de l’étude du mouvement séculaire 4 condi- 
tion que toutes les courtes périodes et notam- 
ment celle d’un an soient correctement filtrées. 
Le mouvement séculaire doit étre calculé a 
partir de toutes les stations ot les variations de 
latitude ont été étudiées assez longtemps que 
pour permettre l’application de la combinaison 
$4S5S6S7. Les variations de longitudes éventuelle- 
ment décelées par les opérations mondiales des 
' longitudes pourraient également étre introduites 
dans ces calculs. 

5. Quelles peuvent étre les causes de ces 
dérives? Constatant l’inexactitude des mouve- 
ments propres adoptés pour certaines étoiles des 
programmes et démontrant surtout que dans 
l’application du programme il existait une dis- 
symétrie systématique marquée entre les di- 
verses stations (fréquence variable.d’observation 
d’étoiles faibles et d’étoiles en fin de'groupe selon 
les observatoires). Melchior avait interprété la 
dérive comme dfie aux inexactitudes des mouve- 
ments propres des étoiles utilisées. Sur la base de 
cette interprétation Phillippov et Sekiguchi ont 
repris indépendamment les calculs de réduction 
et ont montré que la dérive subsistait indépen- 
damment de cet effet. 


DISCUSSION 


Dr. PAUCHENKO* remarked that Dr. Melchior 


said in his communication that until now no one 
has given any explanation of very strong damp- . 


ing of the free nutation of the earth, as found 


* Gravimetrical Observatory of the Academy of Sciences 
of the Ukrainian S.S.R., Poltava, U.S.S.R. 
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Sans doute doit-on considérer que des effe} 
atmosphériques interviennent ici mais d’u1| 
maniére qui reste confuse pour nous (anomalill 
de réfraction liées a la dispostion des group 
d’étoiles par rapport au zénith?). } 

6. Ces restrictions étant faites sur l’exactitud 
des dérives que l’on peut mettre en éviden¢ 
dans le mouvement du péle, on considérera av 
prudence les coincidences a vrai dire troublant} 
que l’on a relevées avec les variations séculair} 
de la vitesse de rotation de la Terre. Sekiguc} 
a remarqué qu’aux deux perturbations les pli 
marquées de la.dérive du pdle moyen: 1922 | 
1936 étaient associées des perturbations aus) 
brusques de la vitesse de rotation de la Te 
Il propose de les attribuer 4 un couplage électr} 
magnétique entre le noyau et |’écorce terrest 

Enfin, faisant abstraction de ces perturbati 
et cherchant 4 déceler seulement l’allure moj} 
enne de la dérive du péle de 1900 4 1950, on co} 
state avec Cecchini que le déplacement du cent} 
de gravité de la polhodie (péle d’inertie Nore 
en 50 ans aurait été de o’091 (2,7 métres) suiv. 
le méridien 37°50’W (direction de Terre Neu 
ce qui serait assez bien expliqué par la fonte 
banquises et des glaciers et la répartition 
eaux sur les océans. 


CONCLUSIONS 


Nous pensons avoir montré ici la complexi 
des problémes qui restent a traiter. Le gain 
précision que l’on obtiendra grace aux nouveal) 
instruments conduira certes 4 de nouveaux pr} 
grés mais seule une trés longue suite d’observ/ 
tions de tous ces phénoménes permettra d 
dégager les lois fondamentales. 
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by Drs. A. M. Walker and Andrew Young. The} 
results may be questioned, not because of tl| 
defects in the method of analysis, but becau 
of the character of the data used. . 
The method of harmonic analysis proposed i} 
Dr. Jeffreys and developed by Drs. Walker ar} 


| 


| 
i 
; fl 1959 April 


r | | Young can be applied only in the case where the 
i'd “coordinates of the pole used for the analysis con- 
. tain no secular parts. The coordinates of the pole 
_used by the International Latitude Service do 
‘not satisfy this condition. The only system of 
‘iy Coordinates which contain only the periodic part 
idl of the movement of the pole is that of Dr. Orlov. 
| Dr. Pauchenko added that he and his col- 
_ leagues had applied a method of calculation 
_ similar to that of Drs. Walker and Young to 
ind analyze the coordinates of the pole in Dr. Orlov’s 
«4 System from 1897 to 1957. They found the time 
», Of relaxation of the free nutation of the earth to 
“be equal to 90 years, and Chandler's period to be 
7 | equal to 432 days. The value thus found for the 
time of liquidation of the free nutation of the 
,. earth provides, apparently, only the lower limit 
, of this magnitude, taking into account that the 
|| ae time and the time of systematic ob- 
servations of the movement of ‘the pole are of 
| the same order of magnitude. 


> mu 
iN 
cent 


a Dr. A. YounG,* speaking of Dr. Melchior’s 


“ second problem, commented on the magnitude 
_of the variation which Dr. Melchior found in the 
Chandler period. Having solved a linear differ- 
ential equation with certain constants (treating 
_the period also as constant), Dr. Melchior found 
_a relation between those constants and the axes 
_of the forced annual elliptical motion. Dr. Young 
li) considered that this was therefore a relation 
il’ between constants and not a functional relation 
\" between variables. He would like to ask if Dr. 
‘1" Melchior had investigated the effect on the 
seri linearity of the differential equation of the vari- 
:/ ation in the Chandler period. Dr. Young said 
_that he had made an, as yet, incomplete investi- 
gation of this using a sinusoidal variation of the 
Chandler period based on Dr. Melchior’s pub- 
lished results and found that there was reason to 
believe that it is not permissible to treat the 
“}) equation as linear if those results are correct. 
If this is so, Dr. Young concluded, then there is 
4 a contradiction in the mathematical treatment 
of the problem. 


veut 
ted 
mn 


8, 


_ Dr. FEpoROov, commenting on Dr. Melchior’s 
Paper, said that he did not wholly agree with the 
Opinion that the apparent variations in the 
Te secular motion of the pole are connected with a 
: * University of Liverpool, Liverpool, England. 
al 
od 
f af 
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change of the International Latitude Service 
observing program. 

If only these observations are taken into ac- 
count, the largest changes in the motion of the 
mean pole are found in 1922 and 1936. It was 
about this time that the JI.L.S. program was 
changed. However, it was also at this time that 
irregularities in the earth’s velocity of rotation 
were detected. Therefore, if one takes into con- 
sideration only the coincidences of time, it would 
be possible to suppose that the velocity of the 
earth’s rotation changed because the I.L.S. 
adopted a new program of observations. 

In addition, Dr. Fedorov thought that at the 
time when the program was changed there were 
in reality no changes in the motion of the mean 
pole (real or spurious) and that Dr. Melchior 
introduced them artificially. Let us suppose that 
we have a continuous curve. Using Dr. Melchior’s 
method we can always obtain a break and, be- 
sides, at any point. If we divide our curve in two 
parts by a straight line x = a, a being an arbitr- 
ary constant, and approximate the left and right 
parts of the curve by straight lines, it is obvious 
that in a general case these two lines will have 
different directions and will intersect the line 
x =a at different points. From this, as we see, 
it does not follow that the initial curve has a 
discontinuity at x = a. 


Dr. MELCHIOR answered Dr. Young by saying 
that the variation of Chandler’s period which 
can be observed has the character of a slow fluctu- 
ation, and that it is not possible to assign a real 
period to it. He added that all methods of 
analysis of the problem gave results in accord- 
ance with those obtained by using the formulas 
given in the second problem. 

To Dr. Fedorov, Dr. Melchior replied that in 
attempting to refute his thesis that the secular 
motions observed for the mean pole are fictitious, 
Dr. Fedorov admitted a priori that the secular 
movement was real. Dr. Melchior noted that the 
transfer of axes of coordinates attained by his 
method coincided with those indicated by 
Kimura and Nicolini by a completely different 
method. Finally, it should be pointed out that 
according to very homogeneous observations at 
Greenwich, there was no drift following the ox 
axis; can one logically admit that drift is found 
only following the oy axis? 
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PART IJ. THE ROTATION OF THE EARTH 
THE IRREGULAR ROTATION OF THE EARTH 


BYE Boule 


Abstract. A mechanical explanation of the observed irregular rotation of the surface of the earth is proposed. From thi 
probabilistic model the correlations between the apparent fluctuations of the motion of the moon in successive years 
derived and a modified least-square method is developed for finding the secular accelerations of sun and moon. 


Let AL be the difference between the observed 
and the tabular longitude of the moon, and 
AL’ the same difference for the sun. Newcomb 
and De Sitter represented AL by a formula 


AE-=x + yl +27? + B, 


where x, y, z are constants, T being the time in 
centuries since 1900, and B an irregular residue. 
Spencer Jones succeeded in representing AL’ by 
a similar formula 


AL! = x! + y'T + 2/T? + B/13.4, 


in which B is the same as before. This means that 
both irregularities may be reduced to zero by 
one and the same clock time correction F?. Thus 
it appears that the irregularities are not due to 
the moon, but to the daily rotation of the earth. 

To explain the irregularities in the earth’s 
rotation it is not necessary to assume (as some 
have done) considerable changes in the moment 
of inertia of the earth. In a paper to be published 
in the Astronomical Journal the following as- 
sumptions were made: 


1. The retarding couple acting on the earth 
due to the friction of the tides is constant = K. 

2. The earth’s moment of inertia is constant 
= C. 

3. The mantle of the earth is a.rigid body. 
Its angular momentum is m = cw, where w is the 
angular velocity and c the moment of inertia. 

4. The earth as a whole is not a rigid body. 
If M is the angular momentum of the whole 
earth, we may define its angular velocity Q by 
ME" CW; 

5. The frictional couple between the mantle 
and the core is proportional to 29 — w; K, = 
f(Q — w). The coefficient f may bea random vari- 
able, varying between a positive lower bound fo 
and o, 

6. In addition to ae there is another couple 
acting on the mantle due to irregular motions 
within the core: Ky = W(T). 

7. This Y(T) is, for any given time T, a ran- 
dom variable with mean value zero; E(W) = 


VAN DER WAERDEN 
University of Zurich, Zurich, Switzerland 


s 


8. Values ¥(7) and W(7>) at.times JT; and f 
far apart are practically uncorrelated. ( 
The differential equations for w and Qare 


“CQ=-K (1 

co = f(Q—w) + W(T) — K. (2 

The solution is ) 
Q =a — pT (3) 
= w) — BT — 0(T), (4) 


where 6(T) is again a random variable with ex: 
pectation zero, and such that values of 6 foi 
times T, and T> far apart are practically uncor 
related. 4 

“Far apart’’ may mean several years or several 
centuries apart; we do not know that beforehane 
and have to find it out from the observations. 
We are not concerned with seasonal variations in 
C, c and w due to melting snow or ice, but only 
with changes in the mean value of w from year to 
year. 

Integrating Eq. (4) from 0 to T, we obtain the 
displacement of any place on the surface of the 
earth (say Greenwich, where most observations 
have been made) since 1900: 


T T . 
ip wdt = wel — 38T? — il 6dt. (5) 
0 


The first two terms give the displacement of a 
fictitious ‘‘mean Greenwich,” which rotates with 
a uniform retardation due to the friction of the 
tides. The last term is the deviation of actual 
Greenwich from mean Greenwich. This devia- 
tion, divided by w, gives the clock time correction 
F*, Multiplying #* by the mean motion of the! 
moon in longitude, we obtain the apparent 


fluctuation of the moon i 
| 


Tr 4 

=k f 6dt (6) 

0 

and the correction to the tabular longitude 
AL =x+ yT + 2T? + B*. 


959 April 
_ The theoretical B* is marked by an asterisk 
o distinguish it from the observed B, which is 
| i by the conventional formula 


a xo + yol + Zl? + Dy, (8) 


F which xo, yo and 2 are fixed values derived 
rom the papers of De Sitter and adopted by 
| aed Jones, Clemence and Brouwer. 

| From Eqs. (6), (7) and (8) we obtain, eliminat- 
ng aE 
| BY = 


{ 


B- (a +f + cI = e [ aar. () 


z 
The problem is to find out whether the con- 
a a, band c can be determined in such a way 
hat 


: *— B-(a+bT +cI) 


x 
es. like an integral k J 6dt, as it should do 


ccording to the theory, to estimate the constants 


i 
| @=x—%, b=y—y, €#=2— 2 


ind to estimate the standard errors of the esti- 
nates of b and c. 

_ Putting ¢ = oin Eq. (9), we obtain a = B(o). 
‘This means that in 1900 mean Greenwich coin- 
‘ided with actual Greenwich. This is an arbitrary 
‘onvention, which may be changed at will. 

_ Subtracting (9) for T + 46T and T — 467, and 


Bh ea ee a 


= Ts 


ontribute to the changes in the earth’s rate of rotation. 


Seven years ago (Brouwer 1952a,b) I at- 
empted a solution of the problem of the irregular 
‘hanges in the earth’s rate of rotation by assum- 
ng that the second differences of the annual 
values of the fluctuation curve in the moon’s 
nean longitude are of a random character, un- 
orrelated from one year to the next. These 
‘econd differences were assumed to have a normal 
listribution with a fixed standard deviation to be 
letermined from the observational data. Their 
nean value would differ from zero by a small 
juantity corresponding to a correction to the 
iecular acceleration term in the moon’s tabular 
nean longitude. 
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dividing by 67, we obtain 


k6 = = — (b+ 2cT). (10) 

Here 6 is the integral mean of 6 over the period 
67, of which T is the middle; 6B is the observed 
change of B during this period. The expectation 
of 6 is zero, hence+ 

b + 2cT is the expectation of the observed differ- 
ence quotient 6B/5T over any period 5T. 

A statistical analysis of the observed 6B’s by 
Van Woerkom shows that differences 6B, taken 
for two one-year periods less than I0 or 12 years 
apart, have a high positive correlation, but if the 
periods are taken 20 or more years apart, the 
correlation practically vanishes. This is just 
what was to be expected from theory. The 
analysis was based upon 130 differences 6B of 
B-values observed in successive years. 

An estimation of b and c from ancient, medieval 
and modern observations by a modified method 
of least squares based upon the statistical prop- 
erties of the 6B’s, resulted in z = 29 +c = 6.1. 

This estimate has a standard error of approxi- 
mately 0.8. The deviations of the observed B’s 
from the parabola B = a + dT + cT? were not 
larger than might be expected from the variance 
of the 6B’s and the errors of the ancient observa- 
tions. 


FLUCTUATIONS AND SECULAR CHANGES IN THE EARTH’S ROTATION 


By DIRK BROUWER 
Yale University Observatory, New Haven, Connecticut 


__ Abstract. Geophysical considerations favor seeking the cause of the random changes in the earth’s rate of rotation in the 
; ‘urbulent motion in the core of the earth. Professor van der Waerden’s analysis applies to this interpretation. The ob- 
ervational evidence indicates that the character of the fluctuation curve may be intermediate between that expected on 
he basis of a theory with frictional couple and one without. If confirmed, this would indicate that two different causes 


The analysis was tested by Van Woerkom 
(1953) in two ways: by examining the properties 
of artificially constructed fluctuation curves and 
by a statistical discussion of the annual values 
of the observed fluctuation curve since 1820. 
He concluded that, on the whole, the observa- 
tional data support the hypothesis, but that 
there may be an unexplained contribution to the 
variance of the first differences. Also, the ampli- 
tude of the observed fluctuations over the 
twenty centuries for which scattered data are 
available appears to be rather smaller than 
might be expected from the character of the 
fluctuations since 1820. Comparison with arti- 
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ficially constructed fluctuation curves shows 
that this situation is improbable but not impossi- 
ble. For the standard deviation of the random 
annual second differences Van Woerkom found 
values between or14 and 022 in time, correspond- 
ing to 0708 to o”12 in seconds of arc in the 
moon’s mean longitude. The uncertainty is pri- 
marily due to the observational errors that 
affect the annual values. 

To some extent these discussions were inde- 
pendent of a formulation of the physical causes 
of the phenomenon. It may be said, however, 
that the assumption of mass movements in the 
earth’s interior, affecting the moment of inertia 
about the axis of rotation, is in agreement with 
the hypothesis. 

Munk and Revelle (1952) presented a compre- 
hensive discussion of the geophysical problems 
that concern us here, and gave strong arguments 
in favor of seeking the cause of the random 
changes in the turbulent motion in the core of 
the earth. Important indirect reasons in favor of 
this explanation are: 


(1) The variations in the westward drift of 
the magnetic field of the earth between 1890 and 
1950 show the proper relationship, qualitatively 
and quantitatively, to the observed changes in 
the length of the day. 

(2) The small displacements of the pole of 
‘rotation are consistent only with nearly sym- 
metrical disturbances, which are difficult to bring 
into agreement with the mass-movement hy- 
pothesis. 


Dr. van der Waerden’s discussion applies to 
the explanation proposed by Munk and Revelle. 
The new element that he introduces is a frictional 
couple proportional to Q — w, the difference be- 
tween the angular velocity of the,earth as a 
whole and the variable angular velocity of the 
mantle. 

The presence of the frictional couple affects 
the character of the fluctuation curve in an im- 
portant way. The first differences of the fluctu- 
ation curve, instead of. being the sum of random 
contributions with mean value zero superposed 
on a linear function of the time, now must be 
satisfied by 


AB=a+6T+ if e(T — t) e-9'dt, 
0 
in which e(T) is a random variable with expecta- 


tion zero. Assume for the purpose of discussion 
that the values of e(T) are uncorrelated from one 
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year to the next. Then the integral may be re 
placed by f | 


€y + e106 * + ene" ae 


Except for corrections at the beginning and th 
end of the interval that become insignificant fc 
large values of n, the sum for 1 years is | 


(E.)/o-en 


The expectation of this sum is zero, but th 
deviations from zero have the mean value «7, 
o is the-r.m.s. value of ¢;. This compares with a 
accumulation of the order on} for the solutio 
without frictional term. Without frictional tert 
the evaluation of the secular acceleration in th 
moon’s mean longitude is affected by a larg 
uncertainty due to the random process. Th 
near-indeterminateness disappears with the it 
troduction of the frictional term. 

The residuals obtained by Van der Waerden 
solution are larger than would be compatible wit 
this assumption. For the moment we can cor 
clude only that with his theory the assumption ¢ 
independent values of « from year to year fai 
to satisfy the observed fluctuation curve. 

However, Van der Waerden’s theory does ne 
require that the values of e are independent fro1 
year to year; he postulates only that «(7) an 
e(T2) at times 7; and T> far apart are practicall 
independent and that the mean value of e over 
long interval of time, which may be sever 
centuries, is zero. 

Pending further study of the problem, it woul 
appear as if the observed fluctuation is inte 
mediate between the two interpretations: of 
with frictional couple, and one without friction 
couple. A not very attractive hypothesis is thi 
the observed phenomenon is the result of tw 
processes acting simultaneously. Munk an 
Revelle remark at the end of their paper the 
organized crustal movement, even if a prié 
unlikely, cannot be altogether ruled out as 
possible cause for the fluctuations in the lengt 
of the day. : 

A more detailed statistical analysis of tt 
fluctuation curve is clearly needed. As I remarke 
earlier, the observational errors in the annu 


values of B since 1820 complicate the preci 


evaluation of the random contributions. TI 


more accurate values that can now be obtaine 


and on which Dr. Markowitz is reporting shoul 
be of great aid in clarifying the issue. 


4} 050 April THE ASTRONOMICAL JOURNAL 99 


i. Pee Fis 
'. Finally, the values of the secular variations in REFERENCES 


Brouwer, D. 1952a, Proc. Nat. Acad. Sci. 38, I. 
. 1952b, A. J. 57, 125. 
Munk, W. and Revelle, R. 1952, M@. N. Geophys. Suppl. 6, 


S33 
Van Woerkom, A. J. J. 1953, A. J. 58, I0. 


VARIATIONS PERIODIQUES ET ALEATOIRES DE LA 
ROTATION DE LA TERRE 


Par NICOLAS STOYKO 
Bureau International de |’Heure, Observatoire de Paris, Paris, France 


_ Résumé. Etude de la periode 1952-1958 a montré la stabilité assez grande de la variation saisonniére de la rotation de 
‘a Terre par rapport a l’heure definitive. Nous avons trouvé 


S = 080240 sin (7 — 31) + 080072 sin 2 (7 + 70), 


yu j est exprimé en jours moyens a partir du 1° janvier. 

| La comparaison avec le temps terrestre des étalons atomiques de fréquence a césium a permis de trouver que la vitesse 
jnoyenne de la rotation de la Terre varie d’une facon aléatoire pendant les trois derniéres années. Pour la période du juin 
‘iu novembre 1955 la durée du jour diminue (la Terre accélere) de 2 X 1075 sec. par mois, du decembre 1955 au mai 1957 
a durée du jour augmente progressivement (la Terre retarde) en moyenne de 4 X 10~* sec. par mois et a partir du mai 
957 elle est presque constante. 

Ainsi, l’extrapolation du temps rotationnel terrestre par rapport au temps atomiques présente des grandes difficultés 


On peut diviser en deux périodes les études des 
rrégularités de la rotation de la Terre de courte 
durée: 


I. Utilisation des horloges a quartz et des 
pendules a partir de 1933. 

_ 2. Introduction de |’étalon atomique a césium 
du National Physical Laboratory 4 Teddington 


| aia par M. L. Essen a partir de juin 1955. 


Quand on utilise les marches des horloges A 
quartz pour l'étude de I'irrégularité de la rotation 
de la Terre, on détermine seulement les termes 
périodiques, dont les périodes ne dépassent pas 
sensiblement une année. Les variations progres- 
sives aléatoires de la rotation de la Terre sont 
mélangées avec les variations progressives de la 
marche des horloges. 

Aprés |’introduction de_l’étalon atomique a 
césium qui est comparé réguliérement avec les 
émissions de Rugby (MSF), il est devenu possi- 
ble de déterminer les variations périodiques et 
aléatoires dans la rotation de la Terre. 

Nous avons repris le calcul de l’irrégularité 
saisonniére de la rotation de la Terre a partir de 
l'année 1952. Les résultats jusqu’au mois de mai 
1955 sont calculés en utilisant les horloges a 
quartz et les heures définitives du Bureau Inter- 
national de l’Heure (B.I.H.). Nous avons utilisé 
les horloges 4 quartz de Greenwich, Neuchatel, 
Paris, Potsdam et Washington. 

Les résultats jusqu’en janvier 1958 sont rap- 
portés a l’heure définitive ou a la moyenne de 
plus de 20 services horaires. Les résultats des 


mois de février et mars 1958 sont rapportés, 
respectivement, 4 la moyenne de I0 et 6 services 
horaires. 

Dans le cas d’étude de l’irrégularité saison- 
niére de la rotation de la Terre en utilisant les 
horloges 4 quartz, il faut déterminer les varia- 
tions progressives de la marche des _ horloges. 
Dans le cas d'utilisation de |’étalon atomique, il 
faut éliminer la variation progressive de sa 
fréquence par rapport au TU2 (Temps Universel 
uniforme provisoire). Cette variation provient 
du changement progressif de la vitesse de la 
rotation de la Terre. Les premiéres études des 
résultats de l’étalon atomique (Stoyko et Stoyko 
1958; Essen, Parry, Markowitz et Hall 1958) 
ont montré que la fréquence de 1|’étalon a césium 
augmente linéairement de 40 Hz par an, c’est a 
dire sa fréquence par rapport 4 TU2 devrait étre 
au mois de mars 1958 égale 4 9, 192, 631,945 Hz, 
étant donné que la fréquence conventionnelle du 
départ (juin 1955) est égale 4 9, 182, 631, 830 Hz. 

L’étude plus détaillée de la variation progres- 
sive de la fréquence de l’étalon a césium a montré 
que sa variation n’est pas linéaire. En prenant 
les différences de fréquence pour I’intervalle de 
12 mois, afin d’éliminer la variation saisonni€re, 
on trouve que ces différences augmentent liné- 
airement jusqu’au milieu de l’année 1957, tandis 
que, aprés cette date, elles diminuent. Cela nous 
montre que la vitesse de la rotation de la Terre 
exprimée en TU2 a des variations aléatoires. 
Dans le Tableau I, nous donnons les variations 
de la durée du jour d’aprés l’étalon a césium 


100 


TABLEAU I. VARIATIONS DE LA DUREE DU JOUR D’APRES 
L’ETALON A CESIUM 


(milieu du mois, en I: 100,000 sec) 


1955 1956 1957 1958 
Janv. +20 + 60 +100 
Féy. +60 +110 +115 
Mars +65 +124 +127 
Avril +48 + 82 
Mai +39 +112 
Juin +06 —I10 + 65 
Juil. —47 —29 + 11 
Aofit —50 <3 + 04 
Sept —2I —05 + 78 
Oct. +20 +74 + 75 
Nov +42 +92 +128 
Déc +46 +84 + 84 


(fréquence conventionnelle 9, 192, 631, 830 Hz) 
par rapport a TU1 (Temps Universel, compte 
tenu du mouvement du pdle). Le signe + indique 
que la durée du jour augmente par rapport a la 
fréquence conventionnelle de l’étalon a césium. 
Nous avons utilisé des méthodes différentes 
pour déterminer la variation progressive de la 
fréquence de |’étalon 4 césium. Dans ce but, 
nous avons utilisé les variations saisonniérés 
déterminées par les méthodes différentes: 


I. extrapolées par le B.I.H.; 

2. déterminées d’aprés les horloges 
pour l’intervalle étudié; 

3. par l’intégration numérique des résultats 
d’aprés |’étalon atomique. 


N 


a quartz 


_ Nous avons fait des combinaisons différentes 
en utilisant ou les valeurs individuelles ou les 
valeurs moyennes pour la période de trois ans. 

On remarque une faible diminution de la 
fréquence de |’étalon 4 césium du juin au no- 
vembre 1955 par rapport 4 TU2 (Tableau II, 
valeur conventionnelle 9, 192, 631, 830 Hz). A 
partir du mois de novembre 1955 la fréquence 
augmente progressivement jusqu’au mois de mai 
1957 avec une vitesse plus ou moins,grande. A 
partir du mois de mai la fréquence reste stable. 

En partant des valeurs de fréquence de |’étalon 
a césium du Tableau II nous trouvons (Tableau 
TABLEAU II. CORRECTIONS DE LA FREQUENCE DE L’ETALON 

A CESIUM 
(milieu du mois) 


1955 1956 1957 1958 
Janv = hw +68 +87 
Fév. +20 +75 +87 
Mars +23 +82 +87 
Avril +27 +85 
Mai +30 +87 
Juin +22 sso tan e7 
Juil. +20 +36 +87 
Aofit +18 +40 +87 
Sept +16 +45 +87 
Oct. +14 +51 +87 
Nov +13 +56 +87 
Déc +14 +62 +87 
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TABLEAU III. VARIATION EROGRESSIVE) DE LA DUREE D 
JOUR D "APRES L’ETALON A CESIUM 


(milieu du mois, en 1/100,000 sec) 


1955 1956 1957 1958 
Janv moo +43 +61 
Féy. —02 +50 +61 
Mars +o1 +56 +61 
Avril +04 +59 
Mai +07 +61 
Juin 00 +10 +61 
Juil. —02 apis +61 
Aofit —04 +17 +61 
Sept —06 +22 +61 
Oct —o8 +27 +61 
Nov. —09 +32 +61 
Déc. —08 +37 +61 


III) les-variations progressives de la durée di 
jour par rapport 4 TU2 du milieu du mois de juin 
On voit que la durée du jour (mai 1957—-matr 
1958) a augmenté par rapport a celle du mois d 
juin 1955 de 0,61 ms. 

En partant des valeurs des Tableaux I et IT] 
on peut trouver les variations saisonniéres de | 
durée du jour, en tenant compte du fait que le 
valeurs du Tableau I sont rapportées a la fré 
quence conventionnelle de |’étalon a césium 
Donc, il faut diminuer les valeurs du Tableau 
de 0,21 ms. 

En faisant la différentiation et l’intégratio 
numérique des valeurs du Tableau I, nous avon 
trouvé deuxiémes valeurs de la variation saisor 
niére de la rotation de la Terre. La moyenne de 
resultats de deux méthodes est donnée dans | 
Tableau IV (d’aprés |’étalon 4 césium par raf 
port 4 TUr; + indique que la Terre retarde). 

Chaque valeur de la variation saisonniére e: 
rapportée toujours 4 la moyenne de 1’année, dor 
la date correspondante est le milieu. 

Dans le Tableau V, nous donnons les résulta’ 
de la variation saisonniére de la rotation de | 
Terre depuis janvier 1952 par rapport 4 TU 
Jusqu’en juin 1955, nous avons utilisé seulemet 
les horloges 4 quartz et aprés cette date | 
moyenne des horloges 4 quartz (poids 1) et d 
l’étalon atomique (poids 2, Tableau IV). 
TABLEAU IV. VARIATIONS SAISONNIERES DE LA DUREE I 

JOUR D’APRES L’ETALON A CESIUM 
(milieu du mois, en 1: 100,000 sec) 


1955 1956 1957 1958 Moy 
Janv +o1 —06 +18 +04 
Fev +41 +41. +33) 5eenso 
Mars +42 -+49,. | -Fa5 anes 
Avril +21 +03 +12 
Mai +09 +30 +19 
Juin —17 —42 —19 —26 
Juil. —68 —63 —76 —69 
Aofit 1 Oy —69 
Sept —=38 —46 —Il —32 
Oct. +09 +23 —03 +10 
Nov +32 +36 +51 +30 
Déc +34 $28 +09 +24 


959 April 


TABLEAU V. VARIATIONS SAISONNIERES 
DE LA DUREE DU JOUR 


(milieu du mois, en 1:100,000 sec) 


1952 1953 1954 1955 1956 1957 1958 Moy 
anv. +27 +or1 +06 +11 +02 —o08 +18 +08 
ee to7e 1031-20) +-3r +38 +36 +33 +28 
Memenerie53 ct 05 1-47 +52 +39 +43. +45 +46 
Kvril ~-+Fi0 +33 +39 +34 +17 +04 +23 
femme t5—4-16 -+13 05 --04 +29 +09 
uin —2I —46 —35 —2I —44 -—20 —31 
ul. —7I —56 -—63 -—66 —60 —70 —64 
loit —37 -—68 -68 —69 -—43 —79 —61 
ete 00) —02)-—30 —34. —4l  —17 —22 
et. +09 +16 +23 +08 +27 +01 +14 
Nov. +26 +19 +18 +34 +38 +45 +30 
Yéc. +10 +12 +22 +27 +28 +17 +21 


La représentation analytique des Tableaux 
‘V et V nous donne pour la variation de la durée 
ju jour (Dj) les formules suivantes, ou 7 est 
ompté en jours moyens 4a partir du Ier janvier 
en ms): 


Ij (52-58) = 0,41 sin (7 + 60,0) 
+ 0,25 sin 2 (j + 115,3) 


Ij (atom) = 0,45 sin (j + 56,7) 
+ 0,24 sin 2 (j + I10,2). 


En partant du Tableau V, nous trouvons les 
orrections dues a l’irrégularité saisonniére de la 
otation de la Terre (Tableau VI). Dans le méme 
ableau nous donnons les corrections d’aprés 
‘horloge atomique (1955-1958, Atom.), les 
orrections extrapolées par le B.I.H. (Exta) 
insi que les corrections extrapolées nouvelles 
Extn). 

La représentation analytique des corrections 
aisonniéres nous donne: 


5 (52-58) = 080240 sin (j — 31,2) 
+ 080072 sin 2 (j + 69,7) 


S$ (atom) = 070259 sin (j — 34,5) 
+ 080071 sin 2 (j + 64,6). 
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Sur les résultats du calcul des corrections dues 
a lirrégularité saisonniére de la rotation de la 
Terre peuvent influencer les erreurs périodiques 
du catalogue d’étoiles fondamentales, les erreurs 
résiduelles dans la variation de marche des hor- 
loges 4 quartz, les erreurs périodiques dans la 
détermination de Il’heure des_ observatoires 
utilisés. 

Dans le cas d’étalon atomique, l’erreur dans la 
variation des marches des horloges 4 quartz est 
remplacée par l’erreur dans la détermination de 
la variation aléatoire de la rotation de la Terre. 
Si l’on admet que la variation de la durée du 
jour a été constante pendant la période 1955- 
1958, l’écart moyen entre les corrections saison- 
niéres individuelles et la moyenne de trois années 
est égal 4 6 ms. En utilisant les variations pro- 
gressives d’aprés le Tableau III, nous trouvons 
pour le méme écart moyen seulement 2 ms. 
Ainsi, on voit l’importance de la méthode utilisée 
pour le calcul des variations saisonniéres. 

Pour apprécier l’influence des différents sys- 
témes sur les valeurs des variations saisonniéres, 
nous avons calculé pour la période de l'utilisation 
de l’étalon atomique les valeurs S dans les sys- 
témes différents (Tableau VII), d’aprés |’étalon 
atomique 1955-1958: 


1. Heure définitive du B.I.H. (Déf) 

2. Heure étalon de 1’U.R.S.S. (URSS) 

3: Heure d’aprés l’observatoire moyen in- 
dépendant du mouvement du pdle (S.P.) 

4. Heure étalon de 1’U.R.S.S. compte tenu des 
corrections dues aux erreurs du catalogue FK3 
d’aprés Nemiro et Tourenko (A. Cat.) 

5. Heure d’aprés l’observatoire de Washington 
(W, deux lunettes P.Z.T.). 


En faisant les différentes combinaisons, deux 
a deux, nous avons trouvé que l’erreur moyenne 
de chaque systéme varie entre 3 et 4 ms. Dans le 


TABLEAU VI. CORRECTIONS DUES A L’IRREGULARITE SAISONNIERE DE LA ROTATION DE LA TERRE 
(debut du mois, en 1:10,000 sec) 


1952 1953 1954 1955 1956 
janv. — 51 — 61 — 80 — 99 — 46 
*év, + 41 — 60 — 63 — 69 — 32 
Mars +146 — 50 — 03 + 23 + 87 
Avril +237 +151 +142 +180 +208 
Mai +257° +253 +266 +284 +259 
juin +202 +303 +310 +299 +269 
fuil. +134 +165 +207 +232 +132 
\ofit — 79 — 5 + 16 + 26 — 51 
sept. —181 —214 —192 —189 —I181 
Jet. —197 —222 — 286 —296 — 304 
Nov. —165 —173 —218 —271 —220 
Déc. — 88 —1I5 —164 —165 —103 
2xta 84 41 19 16 52 
extn 59 33 20 25 31 


1957 1958 Moy Atom Exta Extn 
— 23 — 83 — 63 — 82 —I110 — 70 
— 54 —228 — 38 — 61 — 68 — 16 
steed: + 72 47 Shs) ae oS + 52 
+186 +210 +188 + 198 +160 +160 
+198 +253 + 238 +307 +269 
+290 +279 + 299 +348 +296 
+232 +184 =, 223 +230 +190 
+ 20 =E 12 Se TS ={8 6 = 
—222 —196 = “ey —197 —199 
—275 — 263 — 286 —280 —280 
—272 —220 — 253 —245 —242 
—135 —128 = 158 —160 —149 

36 33 23 sane 24 

28 14 19 24 _— 

26 —— 21 33 14 
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TABLEAU VII. CORRECTIONS DUES A L’IRREGULARITE 
SAISONNIERE DE LA ROTATION DE LA TERRE 


(debut du mois, en I: 10,000 sec) 


Déf URSS S.P. A. Cat. Ww 
Janv. — 82 — 73 — 40 + 04 — 69 
Féy. — 61 —102 + 22 + 10 = 84 
Mars + 58 + 40 +123 +155 + 83 
Avril +198 +190 +214 +263 +245 
Mai +238 +220 +220 +219 +212 
Juin +299 +326 +281 +252 +245 
Juil. +223 +238 +193 +125 +105 
Aofit + 15 + 18 — 33 — 90 + 13 
Sept. —192 —176 —258 —252 —I51 
Oct. — 286 —252 —328 —289 —230 
Nov. —253 —252 —277 —253 —215 
Déc. —158 —180 —118 —144 —160 


Tableau VIII nous donnons la représentation 
analytique des différents systémes du Tableau 
VII, ainsi que de certains autres systémes: 
l’ancienne formule d’extrapolation (Extr anc.), 
Washington pour la période 1952-1958 (Wash 
52-58) et l’heure définitive pour la période 
1952-1958 (Déf 52-58). 

En plus, nous donnons la formule nouvelle qui 
parait représenter mieux que l’ancienne les 
résultats des variations saisonniéres (Extr nouv.). 

On peut tirer les conclusions suivantes du 
Tableau VIII: les amplitudes du terme annuel 
et semi-annuel varient dans les limites assez 
étroites; les erreurs présumées du catalogue FK3 
ne changent pas l’amplitude du terme annuel, 
mais elles changent sa phase de 22 jours environ. 

' Nous n’insistons pas ici sur les causes de vari- 
ations saisonniéres de la rotation de la Terre. 
Les nombreuses hypothéses possibles ont été 
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TABLEAU VIII. COEFFICIENTS DES FORMULES DE 
L’IRREGULARITE SAISONNIERE DE LA 
ROTATION DE LA TERRE 


(en ms) 

Systéme Sint Cos t Sin 2t Cos2t Ama Ams 
Extr anc. +22 —I7 —7 +6 28 9 
Déf. 55-58 +22 =15 =4. SO aera 
URSS 55-58 +21 —-16 -—-5 +6 26 8 
S.P. .55-58 +25 =10 =4 “Geez 7uee 
A. Cat. 55-58 +26 —7 —4 +4 927) 6m 
Wash 55-58. +20) —12) (=—3 [ee beeeeoues 
Wash 52-58 +21 -14 -5 +5 25 7 
Dé. 52-58 +21 —-12 -6 +5 24: 8 
Extr nouv. +22 -13 -5 +6 26 8 |} 


énumerées au Colloque sur les constantes fonda 
mentales del’astronomie en 1950 (Stoyko 1950, 
Les multiples travaux dans ce domaine ont ét 
faits per MM. A. Bilimovitch, J. Cox, F. van de 
Dungen et J. van Mieghem, W. Munk et se 
collaborateurs, N. Pariiskii, N. Stoyko, A. Youns 
Pour le calcul des causes présumées la majorit 
ont adopté les premiéres valeurs de l’amplitud 
du terme annuel de la vitesse de la rotation de | 
Terre. Cette amplitude dépassait largement 2 m: 
A cause de cela ils éprouvaient quelque difficult 
a expliquer l’amplitude de cette variation an 
nuelle. Comme les résultats des derniers calcul 
montrent que l’amplitude du terme annuel n 
dépasse pas I ms, la concordance entre le 
causes et les effets doit devenir meilleure. 
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VARIATION PROGRESSIVE ET VARIATION SAISONNIERE 
DE LA ROTATION DE LA TERRE* 


Par A. DANJON 
Observatoire de Paris, Paris, France 


Résumé. Exposé préliminaire des résultats de la comparaison du temps universel déduit des observations faite 
4 l’astrolabe impersonnel de l’Observatoire de Paris, et du temps uniforme fourni par l’étalon a césium de L. Essen € 


J. V...L. Parry: . 


L’auteur s’est proposé de determiner la vari- 
ation saisonniére de la rotation de la Terre en 
utilisant les observations de temps faites a 
l’Observatoire de Paris avec les astrolabes OPL 1 
ou OPL 15, de 1956,5 a 1958,5, les observations 
étant datées dans l’échelle uniforme de l’étalon 
a césium du Dr. Essen ou dans celle de |’Atomi- 
chron du L.N.R. (Laboratoire National de Radio- 
électricité). 

* Un exposé plus complet et plus détaillé a paru dans les 


Comptes Rendus de |’Académie des Sciences, Paris, Séance 
du 10 Décembre 1958. 


La séparation de la variation saisonniére et d 
la variation lente du temps terrestre est néces 
sairement arbitraire. Il est remarquable qu’o: 
puisse l’effectuer en se donnant a priori les con 
ditions suivantes qui paraissent incompatibles. 

a) On suppose que, dans !’intervalle de 2 an 
considéré, la variation lente peut étre représenté 
par une expression du 3° degré: 


a+ 60+ cé? + d# (6 = t — 1957,00, en années) 


b) On suppose d’autre part que la variatio: 


1959 April 


1957.0 


1957.5 


1958.0 
Figure 1. Variation saisonniére observée, 
1956, 5 4 1958,5. 


Figure 2. Courbe annuelle périodique déduite de deux 
innées d’observation avec la représentation en mnillise- 
ondes 2136 + 936? — 2863 (@ = ¢ — 1957,0) pour la vari- 
ition lente. 


Figure 3. Représentation des écarts. 


saisonniére est périodique, c’est-a-dire que les 
ourbes relatives aux intervalles 1956,5-1957,5 
st 1957,5-1958,5 sont superposables. 

I] était peu vraisemblable qu’une telle décom- 
osition fit possible. Or la solution du probléme 
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a été obtenue sans aucune peine. Les valeurs 
trouvées pour les 3 constantes arbitraires sont, 
en millisecondes: 


6 =-+ 213 


La courbe qu’on obtient alors pour la variation 
saisonniére est assez différente de celle que d’au- 
tres auteurs ont publiée. La différence s’explique 
notamment par l’emploi du catalogue déduit par 
Guinot des mémes observations. 

Les valeurs individuelles étant groupées en 
moyennes, a raison de 2 par mois environ, 
l’écart-type d’une moyenne est: 


¢=+93 d=-— 28. 


C= F2,5 uss 


La représentation graphique des écarts ne fait 
apparaitre aucune marche systématique. 

Ainsi, les résultats de 2 années d’observation 
sont correctement représentés a l’aide de 3 con- 
stantes et d’une courbe annuelle empirique. Ce 
fait est remarquable et il correspond trés proba- 
blement a une réalité physique. 

L’amplitude totale de la variation annuelle 
est de 64 ms. Le maximum de mai est trés aigu. 
peut-étre méme anguleux. I] s’est produit en 
1957,41 et en 1958,41. 

La courbe de la variation saisonniére ne peut 
étre représentée par la somme de deux termes 
sinusoidaux, l’un annuel, l’autre semi-annuel, 
méme en premiére approximation. Les termes de 
période 0,25 et 0,20 sont importants. L’emploi de 
la méthode d’Orlov pour déterminer la variation 
progressive ne parait pas indiqué en pareil cas. 


THE INFLUENCE OF SYSTEMATIC ERRORS OF STAR CATALOGUES ON THE 
DETERMINATION OF THE IRREGULARITIES OF THE EARTH’S ROTATION 


By A. A. NEMIRO anp N. N. PAVLOV 
Pulkovo Observatory, Leningrad, U.S.S.R. 


Abstract. In order to lessen the influence of the errors in the right ascensions of stars of fundamental systems on the 
esults of time determinations, it is proposed that the results of time services observing with transit instruments be used 
or the determination of an independent system of right ascensions of stars. The participation in this work by a maximum 


1umber of observatories of different countries is desirable. 


A detailed study of the irregularities of the 
sarth’s rotation is of great scientific interest and 
yecame possible after the present-day extremely 
uigh precision of measurements of time had been 
ittained with the help of quartz and atomic 
locks. The solution of this problem requires first 
vf all a considerable increase in the precision of 
1stronomical time determinations and a thorough 
exclusion of their errors, especially systematic. 


An important source of such errors are the 
systematic errors of the star catalogue. The 
majority of time services are located in com- 
paratively high latitudes of the Northern hemi- 
sphere and more than half of them between the 
latitudes 50° and 60°. It is just for this zone of 
declinations that the Aa, systematic errors of 
the FK3 are extremely large (Pavlov 1951; 
Nemiro 1958b; Nemiro and Pavlov 1956; 
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Afanasyeva 1957) and therefore noticeably dis- 
tort the system of Universal Time (UT) derived 
on the basis of astronomical observations of time 
services (Nemiro 1958b; Nemiro and Pavlov 
1956; Konstantinov 1955; Smith and Tucker 
1953; Heide 1952). 

In the near future, after the publication of the 

new fundamental system of star positions— 
FK4—this state of conditions will be improved 
considerably. However, it must be noted that in 
spite of the indisputably high quality of the new 
fundamental system—FK4—there is no hope 
that the conditions will be completely satisfac- 
tory. Actually during the past years the precision 
of time measurements with the help of clocks has 
increased much more rapidly than the precision 
of astronomical observations. Especially slow has 
been the increase in the precision of absolute 
catalogues which serve as a base for the compila- 
tion of the reference system of star positions. 
This is connected with the somewhat slow 
progress of the technique of meridian observa- 
tions, the most conservative branch of astrome- 
try. Besides, the compilation of a fundamental 
reference system is very complicated and, to a 
certain extent, a problem open to discussion, 
needing many years for its completion. There- 
fore the difference between the time of publica- 
tion and the mean epoch of observations is 
about twenty to thirty years. 
- It is well known that the star coordinates 
referred to the mean epoch are of maximal pre- 
cision. It follows that when using the fundamen- 
tal system we in fact refer to observations made 
some tens of years ago and naturally they are not 
of the precision which could be attained at 
present. 

The fact that the greatest possible number 
(some tens) of separate catalogues ate used for 
the derivation of the fundamental system does 
not save the situation as many of these catalogues 
are not of high accuracy and have large systema- 
tic errors. The assignment of weights to the 
separate catalogues is an especially controversal 
and complicated problem. Practice shows that 
the precision of fundamental systems, derived 
on the basis of a very limited number of the best 
catalogues, is not inferior in some respects to a 
system derived on the basis of a maximal num- 
ber of catalogues. For example, Eichelberger’s 
Catalogue (for 1925) with respect to the errors 
Aa, and Pu al with regards to the errors Aa, and 
Aue are better than the FK3 in the Northern 
declination zones. 
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Therefore it can be noted that while there is. 
sharp increase in the precision of time measure 
ments with the help of atomic clocks and — 
noticeably increasing precision in the determina 
tion of time with photographic zenith tubes 
Danjon’s prismatic astrolabes and photoelectri 
transit instruments, the precision of fundamente 
catalogues is a hindrance in the solution of th 
problem. 

What is the best possible way out of such 
situation? ( 

The simplest is the smoothing of the righ 
ascensions of the-observed stars on the basis ¢ 
observations made at a time service. This i 
already being done during the reduction c 
observations made with the P.Z.T., where star 
of the magnitude 9 to 10 are used. The righ 
ascensions of these stars are found as a result c 
the smoothing of data of observations with th 
P.Z.T. An analogous process is possible for of 
servations with transit instruments, which ar 
used by most of the time services. An essentiz 
circumstance in this case is that the stars of th 
fundamental system will be smoothed. Therefor 
it will be possible to utilize the observation 
results of time services for improving the funda 
mental reference system. 

This means, first of all, a considerable increas 
in the amount of observational data which wi 
be used for the derivation of the fundamente 
system of right ascensions of stars. If we con 
sider that 77 catalogues were used for the com 
pilation of the FK3 and suppose that all th 
FK3 stars were included in these catalogues an: 
that each star was observed eight times for eac 
catalogue, then we find that each star was ok 
served 626 times. The latter number is probabl 
greatly overestimated. On the other hand, th 
Soviet time services alone, which altogether hav 
18 transit instruments, if observing a commo 
list of stars, could on the average make 2 
observations of each star per year with eac 
instrument. Therefore during two years, 20 > 
18 X 2 = 720 observations of each star coul 
be made and more observational data obtaine 
than those used for the compilations of the FKz 
Besides, the precision of observations wit 
transit instruments of time services can cor 
siderably exceed, with regard to systematic an 


_ accidental errors, those made with large meridia 
instruments. Thus the results of observations ¢ 


time services can be used very effectively fo 
improving the fundamental system of sta 
positions. ) 
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It should be pointed out that the observations 
nade with transit instruments are especially 
uitable for the improvement of the system of 
ight ascensions of the fundamental catalogue 
vith regard to the individual positions of sep- 
irate stars as well as the systematic errors (with 
he exception of the equinoxial error). We would 
ike to point out that the compilation of cata- 
ogues on the basis of photoelectric observations 
yf the Pulkovo Time Service gave very satis- 
actory results (Pavlov 1951; Afanasyeva 1957). 

Observations made with photographic zenith 

ubes are less suitable for improving the funda- 
nental catalogue although they can also give a 
vell smoothed system of right ascensions for the 
tars of the observing list. A complication met 
iere is that all the stars are situated in a very 
1arrow declination zone and besides, these stars 
is a rule are not included in fundamental cata- 
ogues. Therefore their utilization for the im- 
ywrovement of fundamental systems necessitates 
\dditional observations with meridian circles, 
vhich in the given case can hardly give results 
if the highest precision as the azimuth of the 
nstrument must be determined by stars in a 
lifferent declination zone. 
_ The situation is somewhat better if observa- 
ions with prismatic astrolabes are used. In this 
ase, although the process is more complicated 
han that for transit instrument observations, 
he fundamental system can be improved in a 
omparatively wide declination zone, about 50°, 
he precision being the highest in the central 
art of the zone. 

,It must be acknowledged that the transit 
Pe iment i is the most convenient for improving 
he system of right ascensions of the fundamental 
atalogue. It enables an improvement of the 
eference system in a very wide declination zone 
ith about equal precision for different parts of 
he zone. 

_A serious objection to the wide application of 
‘ansit instruments for the time service is the 
slatively low precision of time determinations 
) comparison to that of photographic zenith 
abes and prismatic astrolabes. 

We assume, however, that this objection can 
e€ overcome by a transit instrument of a new 
esign with hermetically sealed pivots, by a 
lore precise determination of the inclination of 
1e axis and the photoelectric registration of the 
loments of star transit. Such an instrument is 
t present being constructed in the shops of the 
ulkovo Observatory. We hope that in the cir- 
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cum-zenith zone it will give the same precision 
of time determinations as the photographic 
zenith tube. A special construction of the pavilion 
is of great importance for increasing the precision 
of the results of observations. Ventilators should 
be used so as to ensure a thorough ventilation. 
Meridian marks are also essential in this case. 
They are especially necessary if the results of 
observations of time services are to be used for 
improving the fundamental system of right 
ascensions of stars, as it is very important for the 
solution of this problem that the observations 
made during one night be of several hours 
duration. 

At present the propositions outlined here can 
be fulfilled by the collective efforts of several 
time services. The greater the number of time 
services taking part in this work, the more rapid 
and more precise the results. International co- 
operation is very desirable. 

In brief, our propositions may be stated as 
follows: that the organization of time deter- 
minations with transit instruments of time serv- 
ices be such that they give a possibility of deriving 
the accidental and systematic errors of the right 
ascensions of the stars of the fundamental cata- 
logue with maximal precision. This can be 
achieved for example if the chain method of 
observations proposed by Nemiro (1958a) 
applied. We point out that this method makes 
possible the exclusion of the systematic errors of 
the catalogue from the results of time observa- 
tions. It is then unnecessary to wait for the 
completion of the final discussion of a large 
number of observations made by time services. 
It is also possible to apply the simpler and less 
laborious method proposed by N. N. Pavlov 
(1958). 

It should be noted that the method of smooth- 
ing star declinations has been used for a long 
time in latitude work. This smoothing gives good 
results because the observations made with 
zenith telescopes are of very high precision. 
This is necessary as the declinations of the faint 
stars used in the Talcott method are influenced 
by large errors. 

Both of these circumstances—the high pre- 
cision of observations and the large coordinate 
errors of the observed stars—are also present in 
time service work and they predetermine the 
expediency of smoothing the right ascensions of 
the stars by observations. 

The use of observations of time services for 
improving the system of right ascensions of the 
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fundamental catalogues proposed here and the 
increase on this basis of the precision of astro- 
nomical time determinations, should, in our 
opinion, be carried out simultaneously with the 
organization at several time services of parallel 
observations with photographic zenith tubes, im- 
personal astrolabes and transit instruments. 
However, we do not deem it expedient that all 
the time services observe only with zenith tubes 
or only with impersonal astrolabes since each of 
these instruments is characterized by its own 
specific errors and since at the same time there 
are still possibilities of increasing the precision 
of transit instruments. 
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Taking into account the above stated, the 
Pulkovo Time Service proposes in the neat 
future to base its work on observations with 
photoelectric transit instruments and a photo 
graphic zenith tube. 
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VARIATIONS IN ROTATION OF THE EARTH, RESULTS OBTAINED WITH THE DUAL 
RATE MOON CAMERA AND PHOTOGRAPHIC ZENITH TUBES 


By WM. MARKOWITZ 
U. S. Naval Observatory, Washington, D. C. 


Abstract. Comparison of photographic zenith tube (P.Z.T.) observations with time derived from quartz-crystal clock 
during 1951 to 1955 and with cesium standards of frequency during 1955 to 1958 indicates that the seasonal variation i 
nearly the same from year to year. Lunar-tidal inequalities of semi- -monthly and monthly periods with amplitudes lo 
about o%001 each were found. A preliminary value of the Love number, k, is derived. Observations made since 1952 wit! 
the dual-rate moon position camera are used to derive AT = ET — ae: Comparison of the P.Z.T. observations an 


atomic standards at the National Physical Laboratory and the Naval Research Laboratory shows details of the iment 


variation from 1955 to 1958. 


1. Introduction. As is well known, there are 
three types of variation in speed of rotation of the 
earth—secular, irregular, and periodic. Although 
the general nature of these variations has been 
known for some time, there is a need for deter- 
mining them in detail. 

Variations in speed are determined by com- 
paring Universal Time (UT), which is based on 
the rotation of the earth, with time based on 
some other standard. This may be’ the moon, 
quartz-crystal clocks, or atomic standards of fre- 
quency. Studies of the variation in speed of rota- 
tion made at the U. S. Naval Observatory have 
been based on observations for UT made with 
the photographic zenith tubes (P.Z.T.’s) at 
Washington, D. C., and Richmond, Florida, and 
for Ephemeris Time (ET) with the dual-rate 
moon position camera at Washington. The clocks 
used have been the quartz-crystal clocks of the 
Naval Observatory and of the National Bureau 
of Standards, and the cesium standards of the 
National Physical Laboratory, Teddington, and 
the Naval Research Laboratory, Washington. 

The problems treated here concern the stabil- 
ity of the seasonal variation, the determination 
of the lunar tidal variation, and the determina- 


‘vantage of the P.Z.T. in the study of periodi 


tion of the irregular variation. The results ob 
tained permit a comparison to be made betwee 
variations in speed determined from astronomica 
observations and those computed from geophysi 
cal and meteorological considerations. A pre 
liminary, but independent, value of the Lov 
number k is obtained, and the nature of thi 
irregular variations is shown. 

2. P.Z.T. The cesium standards are stabl 
to about 1 part in 10. This stability is so higl 
that for all practical purposes the precision wit 
which changes in speed of rotation of the eart! 
can now be obtained depends entirely upon th 
accuracy with which UT can be determined 
The P.Z.T.’s are well suited for the determina 
tion of UT. The instrument is impersonal, th 
zenith is automatically defined by a basin ¢ 
mercury, and no instrumental corrections fe 
azimuth, collimation, or level are require! 
Geen in press). Observations are made a 
the zenith, where refraction anomalies may b 
expected to be a minimum. An important ac 


variations in speed of rotation is that the syster 
of star positions used is made internally cor 
sistent from the P.Z.T. observations themselve; 
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' Thus, periodic errors of star catalogues do not 
‘introduce spurious periodic terms in the deter- 
' mination of UT. 

' The operation of two stations, which are 
| placed about one thousand miles apart, is of 
great advantage, principally because it permits 
_a check on the results to be made. If P.Z.T. No. 
+3, at Washington, and P.Z.T. No. 2, at Rich- 
)mond, both indicate the same change in speed of 
‘rotation a real change may have occurred. If, 
_ however, the changes indicated are not the same 
| the existence of a change is doubtful. 

__ Observations are made on about 200 nights per 
‘year at Washington and 300 nights per year at 
‘Richmond. The Richmond station is advantage- 
/ously located, not only in regard to favorable 
_ weather conditions, but by being near the equa- 
‘tor, at the low latitude: of +25° 37’. At the 
| equator the speed of stars which transit near the 
| zenith is a maximum and the correction for vari- 
‘ation of longitude vanishes. 

| The intercomparison of the P.Z.T.’s shows that 
| systematic errors of the order of 10 milliseconds 
\ may occur, which remain for several days or 
| weeks. However, the cesium standards indicate 
| that the adopted UT2, based on the two P.Z.T.’s, 
is free of noticeable systematic error. For the 
period from 1955.7 to 1958.0 for, example, vy, 
_ the frequency of cesium on UT2, was represented 
“by a straight line. The probable error of a 


‘and that the adopted UT2 was practically free 
’ of systematic error. 
' The suitability of the P.Z.T.’s for determining 
‘ variations of small amplitude is illustrated in the 
‘case of the theoretical lunar tidal terms, whose 
| periods are 13.7 and 27.6 days, respectively. The 
amplitudes of each are about 1 millisecond, and 
‘it was formerly considered that these terms 
‘could not be found from observation. These 
‘terms, however, not only were detected in the 
« P.Z.T. observations of 1951 to 1954, but it be- 
‘came necessary to construct a more detailed 
‘mathematical theory in order to be able to com- 
(pare observation with theory. 
"  P.Z.T. No. 1 was used for time determination 
) at Washington from October 1933 to May 1955. 
i P.Z.T. No. 3, which replaced it at Washington, 
) has been used since May 1954. P.Z.T. No. 2 has 
been used at Richmond since February 1949. 
| The reduction of the Richmond time sights is 
/ made at the Richmond station. The analysis of 
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the results for the two stations, which are tied 
together by monitoring time pulses from WWV, 
is made at Washington. There have been no 
major changes in observing programs or in 
methods of reduction at either station. Hence, 
it is possible to use the P.Z.T.’s to study the 
variations in speed of the rotation of the earth in 
a homogeneous manner. 

3. Seasonal variation. The correction for the 
annual and semiannual terms which comprise 
the seasonal variations may be written, 


ASV =A sin (2at — 6:1) + B sin (4mt — 02) 
= asin 2nt + 6 cos 2zt 
+ c¢sin 4rt + d cos 4rt, 


where A, B, 61, 62, a, b, c, and d are constants, and 
t is the fraction of a year. ASV is to be added to 
UT1, that is, Universal Time corrected for vari- 
ation in longitude, in order to obtain UT2. 

The seasonal variation was first determined 
reliably by Dr. N. Stoyko, at Paris, in 1937. For 
the principal coefficient, A, he obtained about 60 
milliseconds. In subsequent years similar results 
were obtained by Stoyko and by other investiga- 
tors. In 1950 H. F. Finch announced that the 
Greenwich results for 1943 to 1949 confirmed 
those of Stoyko. The seasonal variation thus 
appeared to be stable, with an amplitude of 60 
milliseconds for the annual term. 

In February 1952, however, H. M. Smith of 
Greenwich announced that for 1950 and 1951 A 
had decreased to 33 milliseconds, and for 1952 
Smith and R. H. Tucker obtained 20 millisec- 
onds. The stability of the seasonal variation was 
now questioned. 

The apparent change in the annual term about 
1950 may be due to the fact that quartz-crystal 
clocks previous to this time did not have, in 
general, the reliability of continuous operation 
and precision which was obtained later. In order 
to overcome difficulties due to clock stoppages 
and short runs, various mathematical artifices 
were resorted to, such as working with the mean 
of the second differences of the clock corrections 
of a number of clocks. 

These analytical methods did have the merit 
of showing the existence of the seasonal varia- 
tion, but were probably not sufficiently rigorous 
to furnish definitive values. The analysis cannot 
make up for an intrinsic lack of precision in the 
clocks. 

In 1953 I attempted to determine the seasonal 
variation from the P.Z.T. observations and the 
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clocks at Washington. Let y; be the correction to 
a clock, z. Let y; be that portion of the correction 
due to drift of the crystal, and let o be the sea- 
sonal variation. Then 


Y=Wito. 


It was required that for each clock it should be 
possible to find a function ¥; which would repre- 
sent the drift of the crystal for two years or more, 
and that the o’s obtained from the various clocks 
should be accordant. I was not able in 1953 to 
find any Washington clocks which met these 
requirements. The best clocks were the resonator 
quartz-crystals of the National Bureau of Stand- 
ards, which were placed in operation in May 1951. 
The functions wy; could not be well determined in 
1953 because the initial accelerations were high. 

In October 1954 the resonators were moved to 
Boulder, Colorado. They had then been running 
continuously for more than three years, and it 
was now possible to determine the y;’s. For R2 
and R3 these are, respectively, 


Yo = 31282 — 13.098t1 — 46470¢e-9-1024, (Tt) 


Ws = 18663 — 21.506t; — 18280e-9.001874 


where ¢#; = 0 on August 30.625, 1951. The y; are 
in milliseconds and f; is in days. 

The mean for R2 and R3 of ASV =— a is 
shown in Figure 1. The results of harmonic 
analysis for ASV for individual years are given 
in lines (1) to (3) of Table I. During the 38 
months from August 1951 to October 1954 the 
maximum difference between the o’s for the two 
clocks in accumulated time was never greater 
than 4 milliseconds. 

The Essen ring-crystal oscillator, E-43, which 


TABLE I. ASV, CORRECTION FOR SEASONAL ‘VARIATION 
Unit = o8001 


a b A c d B 
(1) 1950:7-52.7 +21 —25 33 —4 +6 8 
(2). 1952-7-53.7: 84-28) (18 20= 7-8 o. caie 
(3) 1083.7-54-7 th 24 P20, BF Qa ay aie 
(4) 1955 STL Gy Late ot 39 Woe Ole ak a ane 
(5) 1956 7-22"  =55 0 27 = 7 =8 LT 
(6) 1957 Bes aren! ine Stet Svc mir aise EN 
(7) Mean ‘+22.8 —18.0 29.5 —6.8 +7.0 9.7 
(8) Solar tides + 1.5 0.0 —4.3 —1.6 
(9) (7) — (8) +21.3 —18.0 27.8 —2.5 +8.6 9.0 
(10) Winds +17.0 —II.0 20.0 —0.9 +0.4 1.0 
(11) B.LH., 

adopted +22 —17 27.8 —7 +6 9.2 
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Figure 1. The-correction for seasonal variation 
(ASV), 1951.7-1954.7. 


was installed at the Naval Observatory in July 
1954, was used to determine ASV for 1955.0 to 
1956.0. The result is given in line (4) of Table I. 
A preliminary study based on E-56, which was 
installed in October 1954, furnishes a similar 
value of ASV. A more complete analysis will be: 
made later. The drift of E-43 for the interval 
from August 1954 to November 1956, in mill 
seconds, is represented by 


W43 = 56730 — 1.184%. Mi 
— 0.00895f2? + 3.8 X 10 


where ¢2 = 0 on January 1.625, 1956. 

For 1956 and 1957 the solutions were based on 
the running of the cesium standard of the Na- 
tional Physical Laboratory, Teddington (Essen, 
Parry, Markowitz and Hall 1958). The solutions 
are given in lines (5) and (6). 

The mean values of the coefficients are given 
in line (7). Line (11) gives the coefficients which 
were adopted by the Bureau International de 
l’Heure, in 1955, for use during 1956 to convert: 
from UT1 to UT2. The same coefficients were 
retained by the B.J.H. for use during 1957 and 
1958. It is evident that the choice of costicicn ay 
made by Stoyko, was a good one. 

The seasonal variation as determined with the 
P.Z.T.’s is seen to be remarkably consistent from 
year to year with a mean value of 30 milliseconds 
for A. It appears doubtful that there was a sig- 
nificant change in the value of A about the year 
1950, but this possibility cannot be dismissed. 
It would be interesting to know whether formulas 
for drift that are valid for two years or more can 
be found for clocks in use previous to 1950. If 


‘such formulas could be found we might be able 


‘to obtain definitive results for the earlier year. 

4. Causes of seasonal variation. Part of “ 
seasonal variation is due to tides in the crust of 
the earth which are induced by the sun. The 
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‘change in the figure of the earth causes a change 
jin the moment of inertia. Since the angular 
‘momentum is constant the speed of rotation 
changes. Line (8) of Table I gives the theoretical 
corrections for the solar tides Sa and Ssa. Line 
(9) gives the observed ASV when the tidal terms 
are removed, that is, the ASV which is not ac- 
‘counted for by gravitational theory. 

_ The best explanation of the non-gravitational 
annual portion appears to be that it is due to 
‘winds. This explanation was first offered by 
Munk and Miller (1950) in 1950, who obtained 
‘a value of 48 milliseconds for A. In 1951, Mintz 
and Munk (1951), using new data, reduced the 
‘value to 15 milliseconds. Later, in 1953, they 
obtained a slightly larger value of 20 millisec- 
onds, which is not far from the observed value of 
30 milliseconds derived here (1954). The coeff- 
‘cients of the 1953 solution are given in line (10). 
‘In the same year, Pariyski and Berlyand (1953) 
‘also obtained values for A similar to that now 
being obtained by observation. 

_ The observed non-gravitational semiannual 
variation has an amplitude of 9 milliseconds. 
The amplitude of the semiannual effect attributed 
to winds is only 1 millisecond so that at present 
‘there is no adequate explanation of the semi- 
‘annual variation. One might offer a number of 
‘explanations for an annual variation but an 
; unexplained semiannual variation appears to be 
more difficult matter for which to account. 

5. Lunar tidal variations. Tides induced in the 
“crust of the earth by the moon cause variations 
_n speed of rotation in the same manner as the 
sun does. The principal short period terms are 
“lenoted Mf and Mm. The Mf term, due to the 
‘varying declination of the moon, has a period of 
{3.66 days. The Mm term, due to the varying 
listance of the moon ton the earth, has a 
“period of 27.55 days. 

_ The effects of the lunar and solar tidal terms 
/on the variation in speed of rotation were derived 
oy H. Jeffreys in 1928, by F. Andersson in 1937, 
“ind by Mintz and Munk in 1953. The develop- 
“nents were not carried out in detail because it 
was considered that the effects were not capable 
of being observed. 

» After deriving the seasonal variation for the 
yeriod 1951 to 1954 from P.Z.T. observations it 


Che coefficients obtained by observation agreed 
‘vith those computed from theory within the 
‘trors of observation. It was evident, however, 
hat additional theoretical terms were needed. 
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The effective values of the coefficients for any 
lunation depend on the maximum and minimum 
declinations of the moon during the lunation. 

A new calculation of the solar and lunar tidal 
effects was made in 1956 by Woolard (1959). 
The principal terms are as follows: 


Tide Period 
At=k[515.0 sin Q — 18.6 years 
+2.47 sin 2L Mf 13.66 days 
+1.02 sin (2L—Q) 
+2.63 sin g Mm 27.55 days 
—0.17 sin (g+) 
—0.17 sin (g—Q) 
+15.29 sin 20 Ssa 0.5 year 
+ 4.88 sin g’ ], Sa 1.0 year 


where & is the Love number, 0 is ascending node 
of moon, L is mean longitude of moon, g is mean 
anomaly of moon, © is mean longitude of sun, 
g’ is mean anomaly of sun. 

At is the correction in milliseconds to be added 
to Universal Time determined on any night in 
order to remove the lunar and solar tidal effects. 
It should be noted that the lunar terms are 
largely removed by a smoothing process when 
determining the adopted UT2 ordinarily. 

The following process was used to determine 
the lunar terms. Each night of observation gives 
a correction to a clock, E-43 for example, on the 
system UT2. The correction of the clock due to 
drift, 43, is subtracted and the residuals are 
analyzed for terms with arguments 2L and g, 
respectively. 

The coefficients, O, derived from the P.Z.T. 
observations are given in Table II. The values 
for the first 3 years differ only slightly from those 
that were obtained by a different process in 1955. 

The computed coefficients, C, were obtained 
by assigning to k the rounded value 0.300, and 
to @ the value for the middle of the interval 
analyzed. 

From the differences (O — C) an external 
probable error of +0.43 millisecond is obtained 
for an observed coefficient. 

The value of 2 is uncertain. Dr. P. Melchior 
has pointed out that the value obtained from 
earth tide experiments is 0.19 + 0.06 whereas 
the value obtained from the Chandlerian motion 
of the pole is 0.28 + 0.03. There is thus a dis- 
cordance between the terrestrial and astronomi- 
cal determinations. Dr. Melchior noted that 
when enough data have been accumulated we 
may use the lunar tidal terms to obtain an inde- 
pendent value of k. It is interesting to see what 
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TABLE II. COEFFICIENTS OF LUNAR TIDAL TERMS 
Unit = 08001 

(1) (2) (3) (4) 

sin 2L cos 2L sin g cos g 
oO Cc oO (S O (S (0) e 
1951.7-52.7 +1.0 -+1.0 =—0.2 +0.2 +0.4 +0.7 =LOMe 
1952.7-53.7 +1.5 +0.9 +0.5 +0.2 +0.9- +0.7 +0.6 0 
1953-7-54-7 0.7. 0:8 poet i foOnd, +0.5 +0.8 0:20 
1955 =-O:5— 1 On7 0.0 +0.3 +1.1 +0.8 +0.4 0 
1956 +0.3 +0.6 +0.3 +0.3 +2.7. +0.8 +0.5 0 
1957 —0.4 +0.5 +1.0 +0.2 +1.6 +0.9 +0.1 Oo 


value we can obtain with the material at hand. 

We may obtain k from the relation k = 0.300 

O/C. The weight of each determination is C. 
The solutions are: 


k Wt. 

From sin 2L term, 0.270.07, 3.6 
cos 2L term, 0.57+0.21 0.4 

Combined, 2L terms, 0.30+0.07 4.0 
From g term, 0.47 +0.07 357, 
From all terms 0.38+0.05 7.7 


p.e. of unit weight =+ 0.13 


The values of k obtained here are in reasonable 
agreement with the value derived from the 
Chandlerian motion. The value obtained from 
the Mf tide agrees closely, although the Mm 
value disagrees by about twice the probable 
error. It is too early, however, to say whether 
there is a discordance. More data must be ob- 
‘tained, and the effects of water tides on the di- 
- rection of the vertical must be evaluated. 

From a study of the coefficients it is found 
that the observed monthly tide is in phase with 
the tide predicted by theory, with a probable 
error of +1 day. The observed fortnightly tide, 
however, appears to occur 0.5 day in advance of 
the theoretical, with a probable error of +0.5 
day. We may conclude that within the errors of 
observation the observed lunar ‘tides in the 
crust of the earth are in phase with the theore- 
tical. 

6. Irregular variation. The existence of irregu- 
lar variations in the speed of rotation of the 
earth was found, in the past, chiefly from a com- 
parison of the orbital motion of the moon with 
the rotation of the earth. It was not feasible, 
however, to detect changes in speed that occurred 
within a short time in this manner. In fact, it 
was not even certain how these changes took 
place. It was formerly thought that irregular 
changes took place abruptly, but in 1952 
Brouwer (1952) gave another interpretation of 
the irregular. changes. On his hypothesis we 
should expect to observe changes in acceleration 
but not sudden changes in rate. According to 


Brouwer the observed changes in speed may 
the accumulation of small random changes. 
With the development of the cesium standa 
of frequency it becomes possible to determine t 
nature of the irregular variation in detail a1 
choose between the two interpretations. A joi 
program for determining the frequency of cesiu 
in terms of the second of Ephemeris Time is bei: 
carried out by L. Essen and J. V. L. Parry of t 
National Physical Laboratory, Teddington, a1 
W. Markowitz and R. G. Hall of the Nav 
Observatory, Washington (1958). As a prelimi 
ary step, the frequency of cesium in terms 
UT2 was determined for each month from Ju 
1955 to January 1958. The results indicate th 
during this interval the speed of rotation chang 
gradually and not suddenly. From Septemb 
1955 to January 1958 the mean speed of rotatic 
of the earth underwent a practically unifor 
deceleration of 50 parts in 10” per year, equiv 
lent to an increase in the length of the day 
0.43 millisecond per year. There is some indic 
tion of a change in acceleration about August | 
September 1955 but there are not enough data’ 
confirm this. 
7. Moon camera. Ephemeris Time is defin 
by the orbital motion of the earth about the su 
but is obtained in practice from the orbital m 
tion of the moon about the earth. Since 1952 tl 
U.S. Naval Observatory has been using the du 
rate moon position camera to determine Eph 
meris Time (Markowitz 1954, in press). Tl 
camera tracks the moon and surrounding sta 
simultaneously at their respective rates by u 
of a dark plane-parallel glass filter. The image: 
the moon, which passes through this filter, 
shifted in position by tilting the filter during tl 
exposure. It is thus possible to hold the moc 
fixed with relation to the stars. The epoch : 
observation is the instant when the filter 
parallel to a light colored glass plate which is 
. front of the star field. Exposures are for 10 or 2 
seconds. . 
The plates are measured in a machine whit 
has x and y-screws and an accurate rotatir 
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stage. A stroke of a ratchet lever rotates the 
stage 6°. The stars are first measured in x and y. 
The radii of the moon are then measured every 
6° along the bright limb, using only the y-screw. 
The plate is rotated 180° from its original posi- 
tion and the stars are measured again. It is now 
possible to determine the positions of the stars 
and of the center of circle of best fit of the moon, 
with respect to the center of rotation of the stage 
in x and in y. This determines the apparent 
position of the moon with respect to the stars. 

Corrections for parallax are applied to obtain 
the geocentric right ascension and declination of 
the moon. The Improved Lunar Ephemeris, 1952- 
1959, is entered with these quantities and the 
Ephemeris Time is taken out. The Universal 
Time of the epoch of observation is known, so 
that we obtain AT = ET — UT. 

_ The reduction of the moon observation has 
been programmed for the electronic computer 
by my colleague, Dr. R. G. Hall. The reduction 
of all plates taken since the camera was placed 
in operation in June 1952 has been completed. 
The star places used are those taken from the 
Yale Zone Catalogues. Systematic corrections 
N30- Yale, have been determined, but have not 
yet been applied. A provisional correction of 
+0°18 was applied to the AT’s for the difference 
in equinox. 

The results obtained thus far show the exist- 
ence of systematic effects in both the longitude 
and latitude of the moon, which are probably due 
to the systematic departures of the moon from a 
sphere. The period of the libration in longitude 
of the moon is an anomalistic month, or 27.55 
days, and we may expect to find variations in 
AT which are dependent upon the mean anomaly 
of the moon, g. 

The data were divided into periods of 14 
lunations, a period which is nearly-equal to 15 
anomalistic months, and analyzed for corrections 
of the form H sin g + K cos g. The solutions are 
given in Table III. It is probable that the change 
in the coefficients from one period to the next is 
due to the libration in latitude, whose period is 
a nodical month, or 27.21 days. A more detailed 
analysis will be made later, after corrections for 
irregularity of the limb being derived by Dr. 
C. B. Watts have been applied. 


TABLE III. CORRECTIONS FOR TERMS IN g 


Ss 


Lunations Interval H K 

364-381 1952.5-1953.8 +1817 = 0827 
382-395 1953.8 1955.0 —0.08 +0.07 
396-409 1955.0 1956.1 +0.47 —0.06 
410-423 1956.1 1957.2 +0.29 —0.32 
424-437 1957-2 1958.4 —0.23 —0.87 
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TABLE Iv. AT FROM MOON CAMERA 


Epoch Nights ATo ATc 
1952.75 17 29866 29858 
1953-25 16 30.99 30.72 

75 39 30.56 30.32 
1954.25 26 30.22 30.10 
7 24 29.67 29.74 
1955.25 23 30.00 30.22 
of 19 30.24 30.14 
1956.25 23 30.40 30.56 
-75 19 30.26 30.23 
1957.25 23 30.69 30.50 
7 20 31.20 Biss 7 
1958.25 22 31.05 30.96 


Table IV gives the semiannual means of AT. 
AT 0 gives the values initially obtained and AT ¢ 
gives the values when corrected by the solutions 
given in Table III. 

The values of AT¢ are plotted in Figure 2, 
which also shows other solutions for AT. The 
solution by Brouwer (1952), ATz, is a smoothed 
one, based on transit circle observations and oc- 
cultations. The ones marked AT,’ are recent 
values obtained with 6-inch Transit Circle of the 
Naval Observatory. The latest results are in 
press. The ones marked AT,. are from recent 
occultation results. 

The following corrections were applied to place 
the A7’s on a system defined by the equinox of 
N30;ATg (on Newcomb), —1%09; AT,’’ and 
AT. (on FK3), +0%27. 

The curve marked ‘‘AA” gives AA = AT — 
UT, where AT denotes atomic time. This curve 
was obtained by integrating the frequency of 
cesium in terms of Universal Time, previously 
reported (Essen, Parry, Markowitz and Hall 
1958). The two constants of integration, which 
correspond to the value of AT at an initial epoch 
and the frequency of cesium in cycles per second 
of Ephemeris Time, were selected so as to pro- 


Figure 2. Comparisons of UT with Ephemeris and 
Atomic Time (AT). AT = ET — UTand AA = AT — UT 
For symbols, refer to text. 
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vide the best fit with the A7¢’s from the moon 
camera for the interval June 1955 to June 1958. 

A solution of the means of AT¢ from 1954.25 
to 1958.25 gives a value of +0%12 per year’ for 
AT". The cesium standard gives +0816 per year”. 
The agreement is satisfactory. 

The table below gives the values of AT and its 
derivatives from the solution of Brouwer for 
1941.5 to 1950.5, and from the solution based on 
the moon camera and cesium standard for 


1955.5 to 1958.0. 


AT AT’ Ar 
1941.5 23.62 +0843 /year +0%02/year? 
1950.5 28.20 +0.59 +0.02 
1955.5 30.20 +0.20 +0.16 
1958.0 31.10 +0.58 +0.16 


These values indicate that AT” was negative 
at some time between 1950.5 and 1955.5 and that 
two changes of sign occurred. Since it was possi- 
ble to represent the running of clocks R2 and R3 
from 1951.7 to 1954.7 by formulas (1) no very 
drastic changes in acceleration could have oc- 
curred in that interval. The dotted curve in 
Figure 2, marked AT7., represents an estimate of 
the variation in AT from 1950.5 to 1955.5, based 
on all the data now available. A more definitive 
variation will be obtained after corrections for 
the profile of the moon have been applied. 

8. Variations, 1955 to 1958. Since September 
1956 a cesium standard, No. A6, of the Naval 
Research Laboratory, Washington, has also been 
used’ for comparison with UT2 for studying the 
variations in speed of rotation. This standard 
has a quartz-crystal clock associated with it, 
whose frequency is kept matched to that of the 
standard. The cesium beam is not operated con- 
tinuously, but the combination serves as an 
atomic clock. Daily comparisons of astronomical 
time from the P.Z.T.’s versus atomi¢ time are 
made with the Naval Research Laboratory 
whereas only monthly comparisons of frequency 
are made with Teddington. 

The value of the frequency of cesium obtained 
by comparing the cesium standard at the Na- 
tional Physical Laboratory with the observations 
of the moon made at the Naval Observatory 
(Markowitz, Hall, Essen and Parry 1958) is 

ve = 9 192 631 770 cps (of ET). 
Universal Time does not enter in the determina- 
tion of yg, since it cancels out in forming the 
difference 

AT — AA = (ET — UT) — (AT — UT). 


We may therefore use the above value of vz to 
express the length of the day. 
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Figure 3. Variation in length of day, derived fron 
P.Z.T.’s and cesium standards. Ordinate is excess ovel 
86400 seconds of Ephemeris Time in length of day. 
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Figure 3 shows the excess over 86400 seconds 
of ET, obtained by comparing UT1 and UT2, as 
determined with the P.Z.T.’s, with atomic time 
based on cesium standards located at the Na 
tional Physical Laboratory, Teddington and the 
Naval Research Laboratory, Washington. 

The curve marked UT2 gives the mean length 
of the day with periodic terms removed, that is 
the irregular variation. The curve marked UTI 
includes the seasonal and irregular variations 
but not the lunar-tidal terms. 

A rapid decrease of about 0%0003 in the length 
of the day about the beginning of 1958 was indi- 
cated by both the Washington and Richmond 
P.Z.T.’s. There is a possibility, however, that 
what occurred was only a change in acceleration 
about this time. 

It is a. pleasure to acknowledge the help of 
many members of the staff in carrying out these 
studies, in particular, Dr. R. G. Hall. Mr. G. C. 
Whittaker, Director of the Richmond station, 
maintained a continual flow of results. Others 
who participated in making the observations and 
reductions include Miss Ida Ray, Miss Martha 
Sherman, Mrs. Doris Stanley, Mr. J. Siegel, 
Mr. W. Weston, Miss Suzanne Ellis, and Miss 
Norma Grimes at Washington, and Mr. D. 
Monger, Mr. R. Medford, and Mr. J. Rende at 
Richmond. 
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EPHEMERIS TIME 


By G. M. CLEMENCE 
Columbia University, New York, N. Y., and U. S. Naval Observatory, Washington, D. C. 


Abstract. In this article is described the measurement of time from a more general point of view than has been taken 
before. It is shown that mean solar time and ephemeris time are two particular applications of the general principle of 
measuring time by observing angular motions, a principle that can be applied to any angular motion whatever, provided 


that an adequate theory of the motion is available. 


Space and time may be measured in a fashion 
that makes the errors of measurement independ- 
ent. That is, a measure of distance and a measure 
of time may be defined in such a way that an 
error in a measurement of distance will not pro- 
duce any error in the measurement of time, and 
vice versa. It is true that the measurements must 
be restricted to the local frame of reference. 
Special relativity teaches us that the errors are 
no longer strictly independent if two observers in 
relative motion exchange and combine the results 
of their observations. But this restriction does 
not destroy the principle of independence, it 
means only that one observer must have a care 
in interpreting the observations of another. 

It is easily possible to define the fundamental 
units of measurement in such a way that the 
errors in measurements of time and distance are 
not independent. For example, the fundamental 
unit of length might be taken as the meter and 
the fundamental unit of velocity as the velocity 
of light. The unit of time would then be a derived 
one, obtained as the quotient of distance by 
velocity, and then any error in measurements of 
length would be carried over into practical deter- 
minations of time. This procedure would not, 
however, destroy the duality of space and time, 
which is intrinsic and, for any one observer, 
absolute. 

An alternative choice of fundamental units is 
to take the meter as unit of length, and the time 
required for light to travel a distance of one 
meter as the unit of time. In this case also the 
errors would not be independent; practical 
measurements of time-intervals would be affected 
by the error in the length of the meter-bar used, 
as well as by the accidental errors of measurement 
of the time-intervals themselves. 


The non-independence of the errors in the two 
cases mentioned is a result of defining one thing 
in terms of another. Let us consider a third case 
of a different sort. Suppose that the wave length 
of some monochromatic electromagnetic radia- 
tion is taken as the unit of length, and the inverse 
of the corresponding frequency as the unit of 
time. In this case the errors may be independent 
or not, according to the techniques of measure- 
ment employed. If no reference is made to fre- 
quency when measuring lengths, and if no refer- 
ence is made to wave length when measuring 
time-intervals, then the errors are independent. 
On the other hand, the product of a wave length 
by the corresponding frequency is an absolute 
constant, the velocity of light. If the value of 
this constant is known, it then becomes possible 
to infer lengths from measurements of fre- 
quency, and vice versa; and if this technique is 
employed, the errors will not be independent. 
With this choice of units it would become nec- 
essary to examine the techniques used in every 
experiment, in order to ascertain whether the 
errors are independent or not. 

Until now, the fundamental units of length and 
time have always been defined in such a way as 
to preserve the intrinsic independence of the 
errors of measurement. In particular, the recent 
redefinition of the second preserves the inde- 
pendence. It is impossible to infer anything 
about time from measurements of distance, and 
vice versa. I remark, however, that it has been a 
matter of choice rather than necessity. It would 
have been possible to define the unit of distance 
as the distance through which a body would fall 
in a unit of time, under the action of gravity ; and 
if this had been done, the errors of measurement 
would not have been independent. 
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Any angle that is a known continuous function 
of the time, and that can be measured inde- 
pendently of measurements of distance, issuitable 
for the measurement of time. It is not necessary 
that the angle increase uniformly with the time, 
nor even that it always increase with the time; 
it is only necessary that an adequate theory of 
the motion is available. As an example of an 
angle that does not always increase with the time, 
we have the motion of a simple pendulum; but 
in this case we have not a suitable fundamental 
measure of time because we have not an adequate 
theory of the motion, which is unpredictably 
disturbed by variations in gravity, variations in 
the length of the pendulum, and imperfections 
in the suspension. 

The cases where theories of angular motion 
have been brought to the highest state of per- 
fection are the cases of motions of the bodies in 
the solar system. They are of two sorts, motions 
of rotation and motions of revolution. 

Imagine a point fixed on the surface of a rotat- 
ing body, at some distance from the poles of 
rotation. As the body rotates once on its axis, 
this point moves nearly in a circle, describing an 
angle of nearly 2m radians. The equations of 
angular motion of this point are known and may 
be solved, giving 


6 = a+ bt + periodic terms, (1) 


where @ is the amount of angle described since 
some fundamental epoch, ¢ is the time reckoned 
from the fundamental epoch, and a@ and 6 are 
constants of integration. The periodic terms are 
functions of the time and of a and b, and may 
involve additional constants of integration; their 
amplitudes are very much smaller than 27, which 
makes the solution of the equation of motion 
possible by successive approximations. The 
periodic terms express the amounts of precession 
and nutation since the fundamental epoch. The 
number of equations of the type (1) is the same as 
the number of bodies in the system. - 

Next imagine a line drawn from the center of 
mass of a body to the center of mass of the body 
about which it revolves. During one revolution 
of the body this line changes its direction in space, 
describing an ‘angle of nearly 2m radians. The 
equations of orbital motion may be solved for 
this angle, giving 


= A + nt + periodic terms, (2) 


where v is the-amount of angle described since 
the fundamental epoch, ¢ has the same meaning 
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as before, A and 7 are constants of integration 
and the periodic terms now result from th 
gravitational action of other bodies in the system 
The number of equations of the type (2) is on 
less than the number of bodies in the system, th 
sun being excepted. 

The constants, a, 6b, A, and n for each body o 
the system must be determined by observation 
But it is impossible to determine their values fo 
any body until a measure of time has beet 
adopted. In order to define a measure of time it i 
necessary and sufficient to choose numerica 
values for a and_d, or for A and , for one of th 
bodies. Once that has been dene then the re 
maining two constants for that body and th 
four constants for every other body may be deter 
mined by observation. So the number of arbitrarn 
constants of integration for the whole system i 
two less than four times the number of bodies i 
the system (if we neglect the sun and the addi 
tional constants of integration that may b 
present in the periodic terms, and in the expres 
sions for the other components of the complet 
motions). 

In choosing the two constants that define th 
measure of time we are at liberty to take am 
numerical values whatever. The constant 6 or 3 
defines the unit of time and the constant a or £ 
defines the fundamental epoch from which th 
units are counted. 

Consider first the rotation of the earth, and th 
value of b. If we take b equal to 27, then the uni 
of time is approximately one rotation of th 
earth, but not precisely, owing to the presenc 
of the periodic terms in Eq. (1). The precis 
unit of time is one mean rotation of the earth 
the mean value of @ being defined as its valu 
when the periodic terms in Eq. (1) are suppressed 
Again, if we take 6 equal to 47, then the unit o 
time will be two mean rotations of the earth. Or 
if we take b equal to 2 divided by 86400 ant 
multiplied by 1.0027378118868, then the unit o 
time would be the second of mean solar time 
The value of a was determined by the conditio1 
that the mean sun should be over the meridia 
of Greenwich at noon. Having fixed the value 
of a and 0 for the earth it then becomes possibh 
to determine time by observation of the rotatioi 
of the earth, and to determine by observatiot 
the values of a and 6, and of A and n, for othe 


‘celestial bodies. But having adopted mean sola 
‘time as the fundamental measure of time it i 


no longer possible to regard a and 6 for the eartl 
as being constants of integration in the usua 


| 
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sense; their values cannot by any means be 
determined by observations. 

| From what has been said it is now easy to see 
how Eq. (2) for the earth may be used to define a 
measure of time. The values of A and adopted 
were those of Newcomb, because there was no 
reason to prefer any others, and the new defini- 
tion of the second depends on his value of x. 
The values of A and 7 are now not subject to cor- 
rections determined by observations. It is now 
meaningless to speak of an error in the tabular 
mean longitude of the sun, or of an error in its 
tabular mean motion. On the other hand it is 
now possible to determine a and 0 for the earth 
by observations. The measure of time defined 
by Eq. (2) for the earth, with Newcomb’s con- 
stants, is ephemeris time. 

_ Whether Eq. (1) or Eq. (2) is used to define 
the measure of time, the practical determination 


Dr. W. FRICKE,* referring to the paper read 
by Dr. Nemiro, said that the revision of the 
fundamental system of star positions FK3 which 
will result in the formation of FK4 is expected to 
be completed early in the year 1960. At this 
moment a statement is already possible about the 
accuracy of the individual corrections to FK3 
which were published in 1957. According to the 
results obtained from the observations at 
Danjon’s astrolabe at the Observatory of Paris 
during the last few years, these corrections are 
an effective improvement of FK3. The syste- 
matic corrections to FK3, which are being de- 
duced at present by the use of absolute star 
catalogues, will be partly based on catalogues 
‘rom the Pulkovo Observatory, which are ex- 
sellent. Nevertheless, Dr. Nemiro’s remark is 
crue, that also in FK4 the influence of the errors 
of older star catalogues will not be entirely elim- 
nated. Therefore, his proposals have to be 
welcomed that the time determinations with 
nodern transit instruments be organized so that 
they give a possibility of deriving the accidental 
nd systematic errors of the right ascensions of 
he stars of the fundamental catalogue with 
lighest precision. Moreover, the use of the im- 
yersonal astrolabe has to be recommended for 
igh precision observations. 


*Astronomisches echen-Institute in Heidelberg, 


xermany. 
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of time is the same in principle. A table is con- 
structed, showing in adjacent columns values of 
t and of v (or 6). A determination of time consists 
in measuring the value of v (or 8), and consulting 
the table to ascertain the corresponding value of 
te 

It is particularly to be noticed that no assump- 
tion is made about the variability or invariability 
of the measure of time. The only assumption is 
that the equations of motion are known. In the 
case of mean solar time, Eq. (1) for the earth has 
been shown to be inconsistent with Eq. (2) for 
the earth, the moon, Mercury, Venus, and the 
great satellites of Jupiter, whereas the Eqs. (2) 
mentioned are consistent with one another. It is 
this result that has persuaded us that the theory 
of rotation of the earth is incomplete, and to 
discard mean solar time in favor of ephemeris 
time. 


DISCUSSION 


Dr. CLEMENCE thought that in the future the 
Danjon astrolabes should play an important role 
in the construction of fundamental catalogues 
of bright stars. By reason of the numerous ob- 
servations made in the time services these instru- 
ments could greatly reduce systematic errors of 
various sorts in existing catalogues. It would still 
be necessary to rely on meridian circles for the 
equinox and equator, and for the positions of 
reference stars needed in the reduction of photo- 
graphic catalogues. 


Prof. PAVLOV mentioned that in recent 
years the time service is experiencing a great 
change. The function of the International Time 
Service is turning more and more into a service 
for checking the rotation of the earth; at the 
same time, the determination of uniform time is 
becoming in greater degree dependent on the 
atomic standards of frequency. 

In connection with the first problem, the re- 
quirements of the precision of astronomical 
determinations of time, and particularly of 
catalogues of star positions, have significantly 
increased. Prof. Danjon has made a bold attempt 
to deduce seasonal variations of the rotation of 
the earth from observations conducted at a 
single observatory. His work is of great interest. 
The possibility of a similar kind of deduction 
from observations with transit instruments re- 
quires, first of all, a careful elimination of sys- 
tematic errors in star catalogues. 
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PART III. ATOMIC STANDARDS OF FREQUENCY 


ETALONS DE FREQUENCE 


Par B. DECAUX 
Centre National d’Etudes des Télécommunications, Paris, France 


Résumé. De récents progrés ont été obtenus sur les horloges 4 quartz et les étalons atomique de fréquence. 

L’amortissement des cristaux a pu étre notablement diminué par le choix de tailles spéciales et également par l’empl 
de températures voisines du zéro absolu; le vieillissement se trouve en méme temps trés réduit. L’asservissement d’u 
oscillateur par un cristal fonctionnant comme résonateur permet d’obtenir une stabilité accrue. 

Le résonateur 4 jet atomique de césium a fourni de remarquables résultats dans divers pays, et a permis d’affectuer di 
comparaisons internationales montrant un accord a quelques 10~” prés. Des mesures réguliéres sont effectuées par rappo 
aux déterminations horaires. Un de ces appareils est réalisé sous. forme industrielle. 

Le Maser a jet moléculaire d’ammoniac constitue un générateur de courant. I] a donné lieu 4 de nombreuses réalisation 
Certains de ces étalons sont comparés aux horloges d’observatoires, et par l’intermédiaire de signaux horaires, aux étalot 


a césium. Certaines différences ont été constatées. 


Les étalons atomiques de fréquence permettent d’améliorer considérablement la conservation du temps par les horlog 


a quartz. 


Les étalons utilisés actuellement sont basés soit 
sur la résonance mécanique du quartz piézo- 
électrique, soit sur la résonance atomique du 
césium ou moléculaire de l’ammoniac. Les 
cavités résonantes en supraconductivité ne 
paraissent pas étre sorties des laboratoires de 
recherches. 

Pour obtenir la précision la meilleure, il faut 
évidemment que la résonance soit aussi aigtie 
que possible, ce qui se chiffre par le coefficient 
de qualité Q (égal au rapport de la fréquence de 
résonance a la largeur de la courbe pour les 
points d’amplitude égale a 2/2 du maximum). 
Il faut également éviter les influences perturba- 
trices extérieures, dont la plus grave dans le cas 
des quartz est la température, ainsi que les 
phénoménes pouvant faire glisser lentement la 
fréquence dans le temps. 

Diverses améliorations ont été apportées, au 
cours des derniéres années, aux différents ap- 
pareils utilisés, en vue d’obtenir une meilleure 
définition de la résonance et une plus grande con- 
stance de la fréquence. On a pu également réa- 
liser, sous forme industrielle, des étaléns qui jus- 
qu’ici avaient gardé l’allure. d’appareils de 
recherche. 

Rappelons que trois méthodes se présentent 
pour utiliser la résonance, aussi bien d’un quartz 
que d’une raie spectrale: 


1. résonateur. Un générateur auxiliaire est 
réglé manuellement a la résonance, et l’on com- 
pare sa fréquence 4 celle que l’on veut mesurer. 

2. oscillateur asservi. Le générateur auxiliaire 
est réglé automatiquement par un servo-méca- 
nisme, sur la fréquence (ou un sous-multiple) de 
la résonance. 

3. générateur. L’appareil fournit directement 
un courant a la fréquence de résonance. 


QUARTZ 


Les étalons a quartz n’ont pas été, comme o 
pourrait le penser, périmés par les étalons at 
miques. Ces derniers ne constituent en effet, g 
néralement, que des résonateurs, des repéres € 
fréquence, et ne peuvent fournir, tels quels, | 
durée. Il faut donc posséder un garde temps aux 
liaire capable d’assurer la continuité entre k 
mesures effectuées sur les étalons atomique 
Etant donné l’extréme précision de ceux-ci, il € 
nécessaire que les étalons 4 quartz présentent ur 
stabilité aussi élevée que possible, spécialemet 
s’ils ne sont pas asservis au résonateur atomiqu 

Loin de ralentir les recherches sur les quart 
les étalons atomiques ont donc, au contrair 
accentué leur nécessité. Elles ont porté principal 
ment sur trois points: le cristal, la températuy 
la production des oscillations. 

Parmi les tailles du cristal, le type GT | 
l’anneau semblent avoir atteint leur maximu 
de qualité; le barreau et la plaque AT ont < 
contraire fait lobjet de perfectionnemen 
récents. 

Les barreaux, étudiés a l’origine par A. Scheib 
ont été améliorés récemment dans son Institu 
d’abord par une meilleure finition et aussi pi 
l’emploi d’une suspension par ficelage au cent 
d’une ampoule cylindrique portant les électrod 
déposées par métallisation. Les lignes de fore 
du champs électrique étant assez épanouies, on 
disposé, aux extrémités, des électrodes de gard 
réunies a la masse par un condensateur réglabl 
elles permettent l’ajustement fin de la fréquenc 
On a évité l’amortissement qui serait df a uw 
métallisation du cristal lui-méme, ainsi que | 


‘variations dans le temps causées par l’évapor 


tion, la migration des métaux dans la soudur 
etc. Le coefficient de qualité atteint 5a 6 X I 
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a fréquence (100 ou 200/3 kHz) dépend assez 
eu de la densité d’énergie appliquée pour entre- 
enir les oscillations, et le vieillissement reste 
aible. 

Une tendance se dessine actuellement a la 
éalisation de cristaux a fréquence plus élevée, 
le I a 5 MHz par exemple, fonctionnant sur 
artiels, la fréquence fondamentale étant 5 ou 7 
ois plus faible. La taille AT a été ainsi utilisée, 
énéralement sous forme d’un disque circulaire, 
2 plus souvent biseauté pour lui donner un 
spect lenticulaire. Les électrodes peuvent étre 
éposées au centre des faces, ou se trouver sur le 
upport extérieur rigide; la fixation s’opére soit 
ar pingage du bord de la lentille, soit par de 
éritables encastrements. Ces divers procédés 
ournissent a la fois un trés haut coefficient de 
ualité (plus de 10 X 108) et un faible coefficient 
le température. 

Nous citerons principalement dans cet ordre 
7idées les cristaux du National Bureau of Stan- 
ards (N.B.S.), ceux de la Compagnie Marconi, et 
eux qui ont été étudiés en U.R.S.S. Parmi ces 
erniers a été décrite une lentille non métallisée 
tteignant un coefficient de qualité de 17 X 10°. 
Nais une stabilité mécanique meilleure a été 
btenue avec métallisation du cristal et fixation 
ar des soudures situées sur les bords. Un tel 
ristal de 1 MHz, présente aprés polissage un 
oefficient de 9 X 108 dans une atomsphére 
aréfiée d’héhum facilitant les échanges ther- 
iiques. Pour obtenir un faible coefficient de tem- 
érature (aux environs de 70°C), il est nécessaire 
‘assurer l’orientation de la taille 4 1’ prés. La 
Ssistance équivalente n’est que de quelques 
hms et le vieillessement au bout de la premiére 
ia de fonctionnement n’a pas _ dépassé 
xX 10-7. La forme carrée peut aussi étre 
tilisée. 

| Les cristaux classiques. fonctionnent a des 
“mpératures comprises entre +50° et 70°C. 
Jes études entreprises au General Post Office, au 
‘ational Physical Laboratory (N.P.L.) et au 
.B.S. ont montré que la qualité des cristaux se 
Ouvait grandement améliorée pour des tem- 
ratures bien inférieures. En particulier le 
eillissement et |’amortissement se trouvent 
‘duits, déja 4 —10°C. L’emploi d’azote liquide 
’8°K), et surtout l’hélium liquide (2°K), a 
ontré un trés grand intérét; le fonctiounement 
ins I’hélium, actuellement en service régulier 
1 N.B.S., a fourni d’excellents résultats, justi- 
int la complication matérielle de cette méthode. 
‘étude de tailles spécialement adaptées a de 
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telles températures est en cours. On peut aussi 
rappeler les quartz résonateurs du N.B.S. simple- 
ment placés au fond de tubes d’une quinzaine de 
métres forés dans le sol, utilisant la constance de 
la température naturelle. 

L’obtention d’un courant a la fréquence de 
résonance du quartz se heurte 4 quelques diffi- 
cultés lorsque |’on cherche la précision maximum. 
Tout d’abord il faut éviter que les circuits de 
couplage viennent augmenter 1l’amortissement 
que l’on a réussi 4 réduire, a si grands frais. Les 
capacités internes des tubes a vide risquent de se 
modifier avec le temps et avec les tensions 
d’alimentation. Enfin, il faut que l’amplitude du 
cristal reste 4 la valeur minimum, compatible 
avec le niveau de bruit des circuits, pour éviter 
les défauts d’isochronisme et divers autres in- 
convénients. C’est dans ce but que le N.B.S. avait 
utilisé des résonateurs libres; mais ceux-ci ne 
constituaient évidemment que des repéres de 
fréquence. On a donc cherché a asservir, par la 
résonance du cristal, un oscillateur auxiliaire; le 
N.B.S. et la Compagnie Marconi ont réussi a 
obtenir ce résultat sans cristal auxiliaire. La 
Compagnie Marconi vient méme de mettre sur le 
marché un modeéle de ce systéme. Le cristal (a 
5 MHz) est monté dans un pont dont un des bras 
est légérement modifié 75 fois par seconde pour 
faire varier l’équilibre du pont. On recueille ainsi 
a la sortie du pont un signal d’erreur a la méme 
fréquence, dont la comparaison de phase au 
courant produisant la variation du pont com- 
mande un petit moteur; celui-ci entraine un 
condensateur variable de rattrapage connecté a 
l’oscillateur, l’asservissant ainsi, en moyenne, sur 
la fréquence du résonateur. L’étalon Marconi est 
prévu pour rester stable 4 2 X 107” prés, la 
dérive étant compensée automatiquement. 

A. Scheibe et ses collaborateurs ont conservé le 
montage Pierce, dans lequel le cristal est branché 
entre grille et anode d’un tube, en paralléle sur 
un pont de capacités. La stabilité a été trés 
poussée, grace a l’emploi de tensions d’alimenta- 
tion trés réduites (la puissance appliquée a 
l’anode est d’une dizaine de microwatts). 
Monté dans un double thermostat a contacts 
thermométriques et dans un local climatisé, cet 
oscillateur présente un vieillissement de |’ordre 
de 1 X 107°; les irrégularités journaliéres, aprés 
déduction de cette dérive, ne sont que de quelques 
10”, les fluctuations rapides restant dix fois plus 
faibles. 

Des résultats assez analogues sont obtenus en 
France, au Laboratoire National de Radioélectri- 
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cité, avec des barreaux, des thermostats élec- 
troniques a réglage continu, et des oscillateurs 
en pont. Le Laboratoire National de Radioélec- 
tricité (L.N.R.) a également réalisé un oscillateur 
expérimental 4 quartz utilisant exclusivement des 
transistors, y compris pour le fonctionnement 


du thermostat. 


ETALONS ATOMIQUES 


Deux corps sont actuellement utilises dans ces 
étalons: le césium, par sa fréquence de précision 
d’un électron dans le champs magnétique du 
noyau;l’ammoniac, par sa fréquence d’inversion. 
L’un et l’autre ont donné lieu, au cours des 
récentes années, a4 un développement considé- 
rable des recherches et des réalisations. De nom- 
breux appareils sont maintenant en service 
régulier, et comparés aux déterminations ho- 
raires. 

Voici maintenant trois ans que le résonateur 
au césium de L. Essen et J. V. L. Parry est en 
fonctionnement au N.P.L.; il a été décrit en 
détail dans un mémoire de l’an dernier (Essen 
et Parry 1957). Rappelons seulement que les 
atomes issus d’un petit four passent successive- 
ment dans deux entrefers d’aimants puissants, 
qui les trient d’aprés leur moment magnétique. 
Dans l’espace séparant ces deux entrefers, se 
trouvent deux cavités résonantes que traversent 
les atomes, et qui sont excitées en phase par un 
courant a haute fréquence réglé au voisinage de 
la résonance. Un champ magnétique continu 
trés faible, perpendiculaire au jet d’atomes, 
régne entre les cavités. Seuls les atomes en 
résonance atteignent le fil chaud détecteur a 
ionisation, dont le courant décrit la courbe de 
résonance quand on régle le générateur d’excita- 
tion. La fréquence dépendant. du»champ, des 
précautions sont prises pour maintenir celui-ci 
constant et mesurable. Le vide est entretenu par 
des pompes, et l’appareil ne sert que par inter- 
mittence. Des modéles~nouveaux ont, depuis, 
été étudiés. 

La fréquence de transition Fy, (3,0) < (4,0) 
utilisée, a été trouvée par les auteurs égale a 
9.192.631.803 Hz, pour le champ nul, par rapport 
a la seconde de TU2 déterminée par l’Observa- 
toire Royal en juin 1955; N. Stoyko a signalé 
comment cette valeur doit étre corrigée pour 
tenir compte de l’évolution du TU2. Le Dr. 
Essen va nous dire lui-méme les résultats qui 
ont été obtenus et Ja précision que l’on atteint 
avec ce résonateur. 
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La National Company Inc, de Maldc 
(Mass.) construit industriellement 1’‘‘Atom 
chron,” comportant un oscillateur a quartz a 
MHz asservi a la résonance d’un jet atomiqt 
de césium. Le résonateur atomique différe ¢ 
celui du Dr. Essen par quelques particularité 
La plus importante du point de vue pratique e 
la suppression des pompes; le tube (en aci 
imagnétique) est scellé et muni d’un systén 
complexe d’adsorption et d’absorption de gé 
résiduels, maintenant le vide pendant un milli¢ 
d’heures. Les aimants sont extérieurs au tube 
par contre, la_partie centrale de celui-ci form 
guide d’ondés, et alimente en énergie haute fr 
quence les deux cavités résonantes; les cham] 
magnétiques extérieurs sont éliminés par d 
blindages. L’asservissement de l’oscillateur 
quartz s’obtient par un procédé un peu analog 
a celui de l’étalon 4 quartz Marconi et qui ava 
déja été appliqué a la premiére ‘“‘horloge”’ 
ammoniac de H. Lyons. La fréquence d’excit: 
tion du tube a césium, produite par multiplicé 
tion de celle du quartz (en réalité 9.192.631.84 
Hz) est modulé a 100 Hz; le courant ionique a | 
sortie du tube se trouve ainsi modulé, et s 
phase est comparée a celle du courant modulan 
Le signal d’erreur résultant commande un servi 
moteur corrigeant la fréquence du quart 
C’est donc un appareil qui, dans la limite d 
durée du tube a césium, peut constituer u 
générateur a fonctionnement continu; en pratiqu 
il est surtout utilisé pour des comparaisons inte 
mittentes. De nombreux exemplaires de ce 
étalon sont en service aux Etats-Unis, d’autré 
se trouvent actuellement au N.P.L. pour k 
comparaisons dont nous parlera le Dr. Esse 
N. Stoyko et B. Decaux ont réalisé depuis di 
mois la comparaison des étalons du B.I.H. et ¢ 
l’“Atomichron”’ en service au L.N.R. 

Divers facteurs peuvent influencer la precio 
des étalons A césium, en particulier la symétried 
la courbe de résonance, |’énergie introduite dar 
les cavités résonantes, les champs magnétique: 
etc. Les mesures les plus précises exigent donc dé 
précautions relatives a la constance de ces d 
verses influences. 

Ne quittons pas le probléme du césium san 
signaler les recherches effectuées pour exploite 
sa résonance par le procéde du ‘ ‘pompage OF 
tique.’’ De telles études sont poursuivies, sud 


‘autres, par A. Kastler a 1’Ecole Normale Su 
‘rieure de Paris. E 


Depuis quelques années un dévelappemti 
prodigieux s’est manifesté dans l’étude des am 


a 
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olificateurs quantiques, auxquels Townes et 
Sordon ont donné le nom de Maser (Microwave 
Amplification by Stimulated Emission of Radi- 
ition). Sous certaines conditions, des molécules 
xcitées A la résonance par une onde électrique 
‘estituent plus d’énergie qu’elles n’en recoivent, 
éalisant ainsi une amplification extraordinaire- 
nent sélective et exempte de bruit. Comme tous 
es amplificateurs, ceux-ci peuvent, par un mon- 
‘age convenable, engendrer des oscillations qui se 
naintiennent sans excitation extérieure. Parmi 
es divers corps utilisables (dont beaucoup sont 
: ’état solide), l’ammoniac a permis de réaliser 
les générateurs d’oscillations particuliérement 
Peles. Les travaux de Basov et Prokhorov en 
J.R.S.S., de Townes et Gordon aux Etats-Unis, 
le Bonanomi, de Prins, Herrmann et Kartaschoff 
1957) en Suisse, de Shimoda au Japon, ont 
-mené ce Maser a un haut degré de précision. 

_ Il comprend en principe une source d’ammo- 
iac émettant un jet de molécules qui sont 
oncentrées et triées (par une longue lentille 
lectrostatique) ; elles pénétrent dans une cavité 
ésonante de forme allongée, ol elles cédent de 
énergie en s’inversant. Les ondes qui s’entretien- 
ent dans la cavité pour la résonance de la raie 
3,3) ont la fréquence 23.870.127 kHz; leur spectre 
st trés pur et contient trés peu de bruit; leur 
uissance est de l’ordre de 5 X 107” W. 

La lentille de concentration a une forme 
péciale pour éliminer les molécules qui se trou- 
ent dans |’état inférieur d’inversion; elle com- 
rend plusieurs paires d’électrodes cylindriques 
aralléles entre lesquellesest établie une différence 
e potentiel de plusieurs milliers de volts. La 
aleur de cette différence de potentiel réagit, au 
scond ordre, sur la fréquence obtenue et doit 
ts étre maintenue constante; il en est de 
1éme de l’intensité du jet moléculaire. La cause 
1 plus grande de déviation de la fréquence réside 


2 


ans le désaccord de la cavité résonante; il 
itervient habituellement pour 1/1000 dans la 
vabilité générale. La cavité doit donc étre 
‘articuliérement stable, tant au point de vue 
»écanique qu’au point de vue _ thermique. 
€quipe de Neuchatel a réduit trés considérable- 
ent l’inconvénient des déréglages accidentels, 
1 utilisant un systéme de deux cavités réso- 
antes couplées; comme dans le cas des circuits 
'ectriques couplés, le sommet de la courbe de 
sonance de ce systéme présente une partie 
vrizontale plate dans un petit domaine, et par 
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conséquent les variations de réglage de toute 
sorte modifient beaucoup moins la fréquence 
obtenue. Signalons que les mémes chercheurs ont 
réussi a utiliser, au lieu de la résonance de la raie 
(3,3), celle de la raie (3,2), plus faible mais plus 
indépendante des conditions extérieures. 

De trés nombreux Maser ont été réalisés, mais 
il semble que peu d’entre eux aient été utilisés 
pour des comparaisons horaires. Cependant, 
depuis deux ans les horloges de l’Observatoire de 
Neuchatel sont périodiquement mesurées en 
fonction des Maser du Laboratoire Suisse de 
Recherches Horlogéres, et grace aux émissions 
de fréquences étalons HBN la comparaison est 
assurée avec les étalons anglais et francais. 
Citons qu’un modéle industriel de ces appareils 
fonctionne a 1l’Exposition de Bruxelles. D’une 
fagon générale il semble que les Maser, plus 
sensibles aux causes extérieures de variations de 
fréquence, n’aient pas encore atteint un degré de 
stabilité et de reproducibilité aussi élevé que les 
résonateurs a césium. 

En résumé: les divers étalons atomiques de 
fréquence, s’ils ne constituent pas de véritables 
garde-temps, fournissent cependant des repéres 
de fréquence d’une stabilité remarquable, per- 
mettant de suivre, avec une précision instan- 
tanée de l’ordre de 107!° ou mieux, la marche des 
horloges 4 quartz. Baser sur eux seuls une nou- 
velle définition de la seconde serait certainement 
prématuré, mais ils assurent une connaissance 
bien meilleure de l’évolution des garde-temps 
actuels. 
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REPORT ON THE PRECISION OF ATOMIC STANDARDS OF FREQUENCY 


By L. ESSEN 
National Physical Laboratory, Teddington, England 


Abstract. Atomic frequency standards have been checked by laboratory investigations and by the intercomparison 
different standards. 

Existing ammonia Masers enable a frequency to be defined in terms of a spectral line with an accuracy of +1 part 
108 and give a reproducibility under stated conditions of +2 parts in 10. The cesium atomic beam standard enables 
frequency to be defined in terms of a spectral line with a precision of +1 part in 10%. Comparison between the U.K. a: 
U.S. by the radio transmissions of standard frequencies and time signals gave an average agreement of a few parts 


10! but revealed some systematic discrepancies. Direct comparisons between three U.S. standards and the U.K. standa 


gave agreement to 2.2, 3.2, and 1.5 parts in 10%. 


1. Introduction. Any radio frequency spectral 
line is a potential standard of frequency and the 
choice depends largely on practical considera- 
tions. So far two lines only have been used, the 
J = 3, K = 3 line of the inversion spectrum of 
ammonia and the F,m(4,0) @ F,m(3,0) hyper- 
fine splitting line of cesium!. The techniques 
employed in order to approach as closely as possi- 
ble the true line frequency have been described 
and this report is concerned only with the pre- 
cision obtained. The results quoted include some 
work which is in course of publication but the 
information may not be complete as there is 
much further work in progress. 

There are three ways of testing the precision 
of frequency standards: 


(a) The parameters that can vary in practice 
are varied deliberately and the effect is measured 
by comparison with another standard which re- 
mains unaltered. 

(b) Standards of different construction are 
directly compared. 

(c) Standards of different construction are 
compared by means of standard frequency trans- 
missions or time signals. 


All of these methods have been used and it will 
be convenient to discuss the results separately 
for the two types of standard. 

2. Ammonia standards. In the early work on 
ammonia in the U. S. and Japan the absorption 
line was used which was of the order of 10° cycles/ 
second wide at a frequency of 2.3870 X 10° 
cycles/second. Techniques were developed for 
setting on the center of the line and an accuracy 
of 5 X 10-° (100 c/sec) has been claimed ; but it 
is difficult to regard a line having such a width 
as forming a precise standard. A fundamental 
improvement was suggested by N. G. Basov 


and A. M. Prokhorov and by J. P. Gordon, H. J. | 


Zeiger and C. H. Townes, which enabled the line 
width to be reduced to 7 kc/sec.,in the device now 
well known as the Maser. A comparison between 


two similar Masers showed that a stability — 
I X 10-" could be maintained for a short tim 
but it was known from both theory and expe 
ment that the external parameters could pull tl 
frequency from that corresponding to the cent 
of the line. For this reason it has proved difficu 
to define the frequency of operation in terms | 
the spectral line with a precision exceeding + 
part in 108. An extensive investigation has bee 
carried out by J. Bonanomi and his collal 
orators at Neuchatel and they have found that 
a particular Maser is set according to a defini 
procedure a reproducibility of +2 parts in I¢ 
can be achieved. It has been found that some: 
the difficulty arises from the fact that the VJ 
(3,3) line is a multiplet and not a single line. 
Time signals from Neuchatel (HBN) ar 
other signals received at Neuchatel have sin 
July 1957 been expressed in terms of a Maser fr 
quency of 23, 870, 129, 235 c/sec. It has thus bee 
possible to make a comparison with the N.P.1 
cesium standard. The results are given in Table 
The mean difference is of no significance sins 
it depends on the values adopted for the spectr 
line frequencies but the deviations from the mee 
represent changes of the Maser frequency rel 
tive to that of the cesium standard apart fro 
any errors in relating the transmitted frequen¢ 
to the Maser frequency or due to the radio prop: 
gation conditions. | 
3. Cesium standards. The cesium standard % 


TABLE I, COMPARISON BETWEEN NEUCHATEL 
“MASER AND N.P.L. CESIUM STANDARD 


Deviation 

Neuchatel Maser— from mean 

N.P.L. Cs parts value parts 
Date in 10? in 10° 
July 1957 ao SR Paes 
Aug. 1957 = 7-9 +0.4 
Sept. 1957 = G4 +0.9 
Oct. 1957 — 6.4 +1.9 
Nov. 1957 —12.7 —4.3 
» Dec. 1957 — 6.5 +1.8 
Jan. 1958 = zi +0.6 

Mean — 8.3 * 


| 
1959 April 


TABLE II. COMPARISON BETWEEN ATOMICHRONS 
WITH N.P.L. CESIUM STANDARD 


Values of Atomichrons— 
Date N.P.L. Cs parts in ro! 


Atomichron No. 100 (Camden) 


Nov. 1956 — 7 
Dec. —2 
Jan. 1957 — 6 
Feb. — 8 
Mar. — 5 
Apr. Co) 
May —2 
Atomichron No. 106 (Boulder) 
Oct. +7 
Nov. +15 
Dec. +18 
Jan. 1958 +22 
Feb. +15 
Atomichron No. 112 (Harvard) 
Mar. + 2 
Apr. + 1 
May fo) 
June fo) 


the N.P.L. has been in operation since June 1955 
and those constructed by the National Company 
‘Inc., U. S. since the autumn of 1956. A consider- 
able amount of information concerning their 
performance has therefore been accumulated. 
‘The N.P.L. standard was built so that the effect 
of varying the parameters could be studied and 
the results of this investigation have been pub- 
lished. Copies of the paper have been circulated 
to members of Commission 31. It was concluded 
from the investigation that the only two param- 
‘eters that needed careful definition were the 
‘magnetic field and the relative phases of the two 
regions of r.f. magnetic field; and that these 
‘parameters could be readily set to give a pre- 
cision of frequency of a few parts in 10". The 
experimental model served for routine calibration 
of the 100 kc/sec quartz clocks used as working 
tandards with a precision of +1 part in 10 and 
under good conditions the accuracy of setting to 
he cesium resonance is +5 parts in 10". 

The National Co. standard, known as the 
Atomichron, is used with a servo control equip- 
ment to correct a quartz oscillator of nominal 
frequency 5 Mc/sec and maintain it in a fixed 
relationship of 1/1, 838.526368 with the cesium 
frequency. The results of an investigation of its 
performance by the National Co. is in course of 
publication. 

_ It has been compared with the N.P.L. standard 
by means of standard frequency transmissions 
and by time signal comparisons. The former 
ethod can be carried out fairly quickly, meas- 
rements extending over an hour or so giving a 
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precision of +1 part in 10'° but it is necessary to 
take account of variations in the propagation 
conditions. It is believed that over parts of the 
day, the propagation of low frequencies, such as 
the MSF 60 kc/sec, GBR 16 kc/sec and GBZ 
19.6 kc/sec signals is sufficiently stable to give 
an accuracy of +I part in 10! J. A. Pierce 
at Harvard has measured his local standards in 
terms of these transmissions and also in terms of 
a transmission controlled or monitored by an 
Atomichron. In this way the two atomic stand- 
ards have been compared and the results are 
summarized in Table II, in which the average 
monthly values are given. 

The difficulties inherent in such a complicated 
chain of comparison were sometimes severe, and 
the conditions for the comparison with No. 106 
were particularly unsatisfactory the standard 
deviation of these results being about +15 X 
10°. The results with No. 112 which is situated 
at Harvard and thus affords the simplest chain 
of comparison are very satisfactory the standard 
deviation being +3 X 107”. 

Comparisons have also been made by means of 
time signals. In a joint program of work between 
the U. S. Naval Observatory and the N.P.L. 
which was begun in the summer of 1955 in order 
to establish the relationship between the cesium 
frequency and the second of ephemeris time, the 
interval between selected time signals ismeasured 
in terms of cesium at the N.P.L. and astronomi- 
cal time at the U. S. Naval Observatory. Since 
the autumn of 1956 an Atomichron has also 
been measured at the U. S. Naval Observatory 
thus giving a comparison between the Atomi- 
chron and the N.P.L. standard. These results 
are summarized in Table III. 

TABLE III. COMPARISON BETWEEN U.S. NAVAL OBSERVATORY 
ATOMICHRON WITH CESIUM STANDARD 
Atomichron No. 6—N.P.L. 


Date Cs parts in 10 
Sept. 1956 = 6) 
Oct. —I2 
Nov. —I4 
Dec. —18 
Jan. 1957 —10 
Feb. — 7 
Mar. —10 
Apr. — 6 
May oO 
June rae 
July Saat! 
Aug. +10 
Sept. ap te 
Oct. +12 
Nov. + 6 
Dec. +28 
Jan. 1958 +20 
Feb. +11 
Mar. = 4 
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In Tables II and III, an allowance has been 
made for the slightly different magnetic fields 
used in the N.P.L. standard and the Atomichron 
and for the different numbers by which the ces- 
ium frequency is divided. 

They show a somewhat larger difference be- 
tween the two standards than the frequency 
measurements and this discrepancy up to June 
1957 was discussed by the Consultative Com- 
mittee for the Definition of the Second of the 
International Committee of Weights and Meas- 
ures and also at the Boulder (1957) meeting of 
U.R.S.I. It was urged that direct comparisons 
should be instituted in order to eliminate the 
errors possibly associated with radio propagation. 

Dr. Markowitz has also supplied the results 
given in Table IV of comparison in the U. S. 
between different Atomichrons. They probably 
do not explain the discrepancy but they show 
that small differences occur in operation between 
the frequencies of different models. 

During March 1958 a series of direct compari- 
sons were made at the N.P.L., Teddington, 
England, between two Atomichrons, Nos. II! 
and 117, an experimental beam tube 857 XI 
made by the National Co. and the N.P.L. 
standard. Members of the staff of the National 
Co. Inc. and the U. S. Army Signal Engineering 
Laboratories assisted in these measurements. 
The results in Table V give the reproducibility 
of the Atomichron frequencies after resetting the 
electronic circuits associated with the servo mech- 
anisms following a standard alignment proced- 
ure. The N.P.L. quartz clock Q13 was in some 
cases used for interpolation over short time inter- 
vals and the results for this clock given in Table 
V illustrate the general uniformity of drift rate 
obtained. The limits represent the standard de- 


viation of a single observation. 
The Atomichron beam tubes were then tested 
as resonators with the N.P.L. driving source. 


TABLE IV. COMPARISONS BETWEEN 
ATOMICHRONS: IN THE U. S. 


Frequency 
difference 
Date Atomichron designations Parts in 1010 
Nov. 1956 to 6 (NRL)-1 (Camden N. Y.) —3 
June 1957 
Oct. 1957 6 (NRL)-106 (Boulder) +4 
(second half 
of month) 
Nov 6 (NRL)-106 (Boulder) Co) 
Dec. 6 (NRL)-106 (Boulder) fe) 
Jan. 1958 6 (NRL)-106 (Boulder) +3 
Feb. 6 (NRL)-106 (Boulder) -3 
Mar 6 (NRL)-106 (Boulder) —4 
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TABLE V. FREQUENCY REPRODUCIBILITY 
AFTER .REALIGNMENT 


Unit-1 <c1opx? 


Frequency difference 
(857 X 


+1II see Q13- 
Date 1958 rr1-N.P.L. 117-N.P.L. N.P.L. N.P.L. 
Mar. 6 3.6 3.5 230.0 
7h 3.8 3.6 231.8 
10 3.3 233.5 
II 4.0 233.6 
12 3.2 4.7 234.5 
14 2.5 4.2 235.3 
17 4.1 3.5 238.2 
18 (a.m.) 3.7 3:7 239.4 
18 (p.m.) 3.6 27 238.1 
19 (a.m.) 3 ang 240.2 
19-(p.m:) 4.0 3:2 240.0 
20; (a.m:)/= 13:9 3:1 240.7 
20 (p.m.) 4.5 241.4 
24 (a.m.) 3.5 3-4 244.7 
24 (p.m.) 3.6 2.8 244.8 
24 (p.m.) 3.8 3-3 245.5 
26 3.4 244.1 
28 a 245.8 
Apr. 1(a.m.) 3.5 3.0 249.2 
1 (peml) eee7) 2.9 249.6 
3.5 “2 250.7 
Mean 3.740.3 3.2+0.4 4.50.3 


Owing to the greater length of the Atomichro 
tubes and their resulting lower band widths the 
could be used with slightly greater precisio 
than the N.P.L. beam tube and it is believe 
that a setting accuracy of +2 parts in 10" we 
achieved. The measurements were used in a 
attempt to analyze the effects that could caus 
the frequency to depart from the true line fre 
quency, ‘and these results are reproduced i 
Table VI. 

The measured errors due to the servo, elec 
tronics and asymmetry can legitimately be af 
plied as corrections to the mean frequencies « 
Table V and the final results of the comparison 
are then as follows: 


11I-N.P.L. 2:2 1.4K 
117-N.P.L.. 332) 1:4, Om 
857 XI-N.P:L. 1.5 22 1.4 X Dor? 


4. New developments. This subject is develo 
ing so rapidly that a new and important advanc 
has been disclosed since the first draft was pre 
pared. If a buffer gas, such as argon at pressure 
up to 30 mm Hg is included with cesium vapoi 
collisions of the cesium atoms with the buffer ga 


TABLE VI. BEAM TUBES TESTED AS RESONATORS 
Unit 1 X 1074 , 


Frequency error : 
857 XI 
+11Ie N.P.L. 
Servo & electronics +8+4 + 6+4 +20+4 O:= 
Asymmetry —3+7 —16+7 o+7 —I10+4]1 
R.f. power +4+2 + 942 O+2 ot 


Effect Ilr 117 


959 April 


1o not produce transitions: but greatly reduce 
he average velocity in a particular direction, and 
1ence the broadening of a line due to the Doppler 
ffect. An optical means of detecting transitions 
yetween the hyperfine levels has given an in- 
reased sensitivity and a line 40 c/sec wide (in 
),192 Mc/sec) has been observed. The actual 
requency is influenced by the buffer gas, however, 
ind the new techniques do not at present pro- 
ride a definitive standard. 

5. Conclusions. It is concluded from this review 
hat the Maser-type ammonia standard can be 
ised to define frequency in terms of a spectral 
ine with an accuracy of +I part in 108 and that 
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the reproducibility of a carefully defined Maser 
is considerably better. 

Existing cesium standards can be used to 
define frequency in terms of a spectral line with 
an accuracy of +2 parts in 10! and the results 
have been integrated throughout a time interval 
of three years for the purpose of comparing 
atomic and astronomical times. 

There are goods reasons for believing that the 
accuracy of the atomic beam-type of cesium 
standard can be increased by at least ten times 
and that experiments with new techniques such 
as optical pumping, and the use of buffer gases 
will lead to simpler forms of standard. 


DISCUSSION 


Dr. STOYKO commented that even though 
he atomic standard is not a clock, it can still be 
ised as a time-keeper through the intermediary 
ff quartz clocks. It has been possible, through 
he kindness of Dr. Decaux, to use Atomichron 
No. 107 at the National Laboratory of Radio- 
lectricity in Bagneux to interpolate from the 
itomic time at the Time Service of the Observa- 
ory of Paris and the International Bureau of 
fime. This atomic standard is compared each 
lay with three quartz clocks which, in turn, are 
‘ompared with the Time Service at the Observa- 
ory of Paris. The integration of frequency of the 
ltomic standard in conjunction with these three 
juartz clocks permits control of the precision of 
ime maintenance by this method. Since March 
958 this method has been used for the direct 
aterpolation of the correction of the main clock 
«cording to the atomic observations of the 
Ibservatory of Paris. Since the conventional fre- 
juency of the atomic standard does not corre- 
pond to the exact second with the provisional 
Jniform Universal Time (UT2) which is the 
yasis of time services, the correction of the atomic 
lock must be represented as calculated by a 
‘near function which is of a sufficient precision 
r an interval of the order of two months. 


Prof. MIKHAILOV, the Chairman, asked 
ow long a cesium standard would run. 


Dr. ESSEN reported that in the United 
ingdom it has been arranged to use the stand- 
rd which controls the standard frequency trans- 
issions for generating the signals for the time 
rvice of the Royal Greenwich Observatory. 
f the frequency remains as at present based on a 
alue of 9,192,631,830 c/sec for the cesium reson- 


ances step adjustments of about one per month 
would be made to the signals to keep them rea- 
sonably close to UT2. If the frequency were 
changed to 9,192,631,770 c/sec, the value in 
terms of the second of ephemeris time still more 
frequent adjustments would be required. It is 
therefore proposed to operate the service in the 
following way: 


The frequency of 9,192,631,770 c/sec will be 
used as the standard but the service will operate 
with an error of about +160 parts in 10! so as 
to bring the frequency and time pulses in close 
accord with the time scale of UT2. Corrections 
published by the Royal Greenwich Observatory 
will give the errors from UT2. Corrections pub- 
lished by the National Physical Laboratory will 
give the errors from the cesium frequency of 
9,192,631,770 c/sec to +1 part in 10%. This will 
also be the error from ephemeris time with the 
lower precision of +20 parts in 10". All possible 
time and frequency information will thus be 
given by the service. 


Dr. MARKOWITZ commented that the 
operation of a cesium clock does not require that 
the cesium beam itself shall be in operation con- 
tinuously. The cesium clock at the Naval Re- 
search Laboratory, Washington, D. C., consists 
of a cesium beam frequency standard which is 
operated in conjunction with a quartz-crystal 
clock of high precision which indicates time. 
The quartz-crystal clock is operated continu- 
ously, and is adjusted in frequency daily to that 
of the cesium standard. It is not necessary, 
therefore, to operate the cesium beam continu- 
ously, and in fact, the beam itself may be 
changed. 


(End of Symposium) 
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